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PREFACE 


This  report  presents  an  analysis  of  the  libration  and  vibration  of  a 
passive  communications  satellite  known  as  PACSAT,  Previous  related  Rand 
reports  include  N-1780-ARPA, ^PACSAT:  ~A~  Passive  Communication  Satellite 
for  Survivable  Command  and  Control,  November  1981,  and  R-2920-ARPA, 
Electromagnetic  Properties  and  Communicat ion  Characteristics  of  PACSAT  ; 
August  1982. 

The  research  is  sponsored  by  the  Strategic  Technology  Office, 

Defense  Advanced  Research  Projects  Agency,  under  the  supervision  of  Dr. 
Sherman  Karp.  It  is  an  element  of  a  broad  range  of  strategic 
communications  studies  being  conducted  for  ARPA  by  Rand. 
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SUMMARY 


The  passive  communications  satellite  known  as  PACSAT  has  been  under 
consideration  for  many  years.  The  satellite  could  be  useful  for  command 
and  control  of  strategic  forces  because  it  can  provide  survivable 
communications.  Many  studies  have  shown  that  it  can  provide  low-to* 
medium  data  rates  at  distances  up  to  a  few  thousand  miles  between 
moderate  sized  mobile  terminals. 

PACSAT  consists  of  a  uniform  linear  array  of  electromagnetically 
resonant  spheres  which  provides  a  very  large  scattering  cross  section  in 
a  particular  direction.  It  is  nominally  vertical  and  straight,  and  is 
maintained  in  that  condition  by  gravity-gradient  stabilization.  Previous 
studies  have  shown  that  verticality  must  be  maintained  to  1.5  deg,  and 
straightness  to  1/8  wavelength  (.5  cm  at  the  desired  operating 
frequency) .  This  report  addresses  the  effects  of  various  disturbing 
mechanical  forces  on  verticality  and  straightness.  Certain  of  these 
effects  are  quite  strong,  and  it  Is  shown  in  particular  that 
micrometeoroid  impacts  and  thermal  bending  effects  cause  motions  which 
far  exceed  the  permissible  displacements. 

The  general  problem  is  formulated,  with  gravity-gradient  tension  and 
elastic  forces  regarded  as  the  principal  forces,  and  other  forces  treated 
as  perturbations.  The  motion  is  separated  into  center-of-mass  motion, 
rigid-body  libratlon  around  the  local  vertical,  and  vibration  transverse 
to  the  libration.  The  fundamental  modes  of  PACSAT  in  a  circular  orbit 
with  no  external  forces  are  studied  in  detail.  The  motion  is  stable, 
undamped,  and  maintains  the  original  amplitude,  so  the  verticality  and 
straightness  requirements  must  be  satisfied  by  the  initial  conditions,  a 
very  difficult  task. 

The  effect  of  orbital  ellipticity  is  the  first  perturbation  studied. 
It  Induces  an  in-plane  libratlon  of  amplitude  in  radians  equal  to  the 
eccentricity.  The  out-of-plane  libratlon  is  unstable,  but  the  exponent 
is  fourth  order  in  the  eccentricity,  and  the  instability  takes  thousands 
of  years  to  develop.  There  is  no  significant  effect  on  vibration. 


Earth  oblateness  effects  are  studied  next.  The  precession  of  the 
orbital  plane  couples  the  in-plane  and  out-of-plane  motions,  producing  a 
resonant  excitation  of  the  out-of-plane  motion,  with  an  amplitude 
proportional  to  the  product  of  eccentricity  and  inclination  for  near¬ 
circular,  near-equatorial  orbits.  Despite  the  resonance,  the 
displacement  is  small.  There  are  no  discernible  effects  on  vibration. 
Lunar  and  solar  gravity  and  solar  radiation  pressure  all  prove  to  be 
negligible. 

Micrometeoroid  Impacts  produce  important  vibrations.  ft  is  shown 
that  a  particle  of  a  few  micrograms  mass  at  meteoric  velocity  can  cause 
deflections  which  exceed  the  straightness  requirements,  and  the  impact  of 
such  a  particle  can  be  expected  as  often  as  once  a  month.  Even  if  the 
satellite  is  initially  vertical  and  straight,  it  will  not  be  able  to 
maintain  Its  condition  under  expected  meteoric  bombardment. 

Thermal  bending  proves  to  be  the  most  important  cause  of  vibration. 
The  solar  heat  input  causes  differential  expansion,  which  makes  PACSAT 
curve  away  from  the  sun.  Eclipsing  by  the  earth,  or  shadowing  within 
PACSAT,  causes  the  solar  input  to  vary  so  as  to  excite  the  lowest 
vibrational  frequency  of  PACSAT,  which  is  at  approximately  three  times 
orbital  frequency.  The  resonance  causes  large  displacements.  Unless  the 
absorptivity  of  the  surfaces  is  kept  below  .0007,  which  is  well  beyond 
the  current  state  of  the  art,  the  straightness  requirements  cannot  be 
met . 

It  is  concluded  that  the  flexural  misbehavior  of  PACSAT  in  orbit  is 
such  that  it  is  most  improbable  that  the  present  design  {unsupported 
linear  array)  can  perform  its  communications  function. 
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I.  INTRODUCTION 


The  survivability  and  endurance  of  the  coraraunications  systems  used 
for  the  command  and  control  of  strategic  forces  are  vital  issues. 

Concern  that  existing  systems  may  be  vulnerable,  despite  their  diversity 
and  specialized  design,  has  led  to  numerous  suggestions  for  their 
improvement.  Among  these  is  a  proposal  by  Yater[l]  for  an  advanced 
passive  communications  satellite  concept.  If  it  proves  practical,  this 
design,  now  called  PACSAT,  may  be  useful  for  the  command  and  control  of 
strategic  forces  because  it  can  provide  survivable  point-to-point 
communication  at  low-to-medium  data  rates  at  distances  up  to  a  few 
thousand  miles  between  moderate  sized  mobile  terminals. 

The  PACSAT  concept  was  investigated  at  the  Stanford  Research 
Institute[2] ,  and  in  a  recent  treatment  by  BedrosianlS] .  PACSAT  consists 
of  a  uniform  linear  array  of  scatterers  whose  axis  is  nominally  pointed 
toward  the  center  of  the  earth.  Gravity-gradient  stabilization  is  used 
to  keep  the  array  erect.  When  illuminated  by  a  transmitter  on  the  earth, 
the  array  scatters  energy  much  like  a  diffraction  grating.  If  the  angle 
of  incidence  (the  angle  between  the  line  of  sight  from  the  transmitter  to 
the  array  and  the  array  axis)  is  denoted  by  a,  the  first-order  grating 
lobe  is  directed  back  toward  the  earth  at  a  half -angle,  B,  given  by  the 
relation: 


s(cosa  +  cosP)  =  X 


1.1 


where  s  is  the  element  spacing  and  X  is  the  wavelength. 

The  region  illuminated  by  the  grating  lobe  on  the  earth  will  be  a 
narrow  annulus.  A  communication  link  is  established  by  choosing  the 
wavelength  so  the  annulus  includes  the  desired  receiver.  The  radar  cross 
section,  available  bandwidth,  and  other  communication  parameters  are 
calculated  in  Ref.  3,  where  it  is  also  shown  that  coverage  considerations 
lead  to  placing  PACSAT  in  a  nominally  circular  equatorial  orbit  at  an 
altitude  of  1-1.5  earth  radii. 
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Reference  3  calculates  the  limits  on  PACSAT  performance  caused  by 
departure  from  the  ideal  configuration.  Effects  considered  are  near¬ 
field  effects  (PACSAT,  at  a  length  of  1,500  m,  a  wavelength  of  3.75  cm, 
and  the  indicated  altitude,  is  in  the  transition  zone  between  Fresnel  and 
Fraunhofer  diffraction),  librational  motion  about  the  vertical,  and 
distortion  of  the  array  from  perfect  straightness  into  a  parabolic  shape. 
The  reduction  of  the  radar  cross  section  is  calculated.  Upon  extracting 
the  information  from  the  extensive  curves  of  Ref.  3,  the  requirements  for 
vertlcality  and  straightness  are  obtained  and  shown  in  Table  1. 

Thus,  if  the  radar  cross-section  loss  is  not  to  exceed  1  dB,  the 
array  must  maintain  verticality  to  within  1.5  deg,  and  must  be  straight 
to  about  0.5  cm  at  the  wavelength  of  interest.  For  an  array  length  of 
1,500  m,  if  the  parabola  be  approximated  by  an  arc  of  a  circle,  the 
radius  of  curvature  of  the  circle  must  exceed  50,000  km,  about  eight 
times  the  radius  of  the  earth.  This  is  clearly  a  very  difficult  task, 
since  the  straightness  must  be  maintained  in  the  presence  of  a  variety  of 
disturbing  forces.  If  the  distorted  shape  were  other  than  parabolic,  the 
limitation  would  be  applied  to  the  rms  displacement.  The  purpose  of  this 
report  is  to  investigate  the  phenomena  and  ascertain  whether  PACSAT  can 
meet  the  flexural  requirements . 

The  behavior  of  PACSAT  in  orbit  is  closely  related  to  the  theory  of 
the  attitude  dynamics  of  satellites  with  flexible  appendages.  This 
subject  has  received  extensive  treatment  in  the  literature.  A  review 
article  by  Modi [4] ,  written  in  1974,  lists  over  200  references,  and  there 
have  been  numerous  publications  since.  Although  PACSAT  at  1,500  m  length 

Table  1 

PACSAT  FLEXURAL  REQUIREMENTS 


Reduction  in  Departure  from  Maximum  Parabolic 

Cross  Section  (dB)  Vertical  (deg)  Displacement  (wavelengths) 


1 

3 

10 


1.5 

2.0 

2.7 


.  127 
.217 
.401 


is  longer  than  any  satellite  ever  flown,  it  is  most  similar  to  the  Radio 
Astronomy  Explorer  satellite  (1968-55A),  which  was  placed  in  a  5,850  km 
altitude  circular  orbit  in  July  1968  and  carried  four  228  m  extensible 
booms.  The  dynamics  of  that  satellite  have  been  well  studied[5,  6]. 

More  recently,  the  motion  of  a  gravity-gradient  stabilized  satellite 
driven  by  orbital  eccentricity  has  been  treated  by  Hablani  and 
Shrivastava[7]  (see  also  comments  by  Alfriend[8] ) ,  and  the  stability 
diagram  for  the  Mathieu  equation  which  governs  the  dynamics  of  a  flexible 
body  has  been  Investigated  by  Kumar  and  Bainum[9]  (see  also  comments  by 
MacNeal(10]). 

The  dynamics  of  PACSAT  in  orbit  involve  the  motion  of  the  body 
arising  from  at  least  the  following  causes: 

1.  Earth's  central  gravitational  field 

2.  Elastic  forces  of  tension  and  bending 

3.  Effects  of  initial  conditions 

4.  Ellipticity  of  orbit 

5.  Earth  oblateness  effects 

6.  Solar  and  lunar  gravity 

7.  Solar  radiation  pressure 

8.  Micrometeoroid  impacts 

9 .  Thermal  bending 

The  first  two  of  these  will  be  regarded  as  fundamental,  and  all  of 
the  rest  will  be  treated  as  perturbations.  It  will  be  shown  that  the 
most  important  effects  which  cause  libration  about  the  vertical  are 
initial  conditions,  ellipticity,  and  oblateness,  while  the  roost  important 
causes  of  departure  from  straightness  are  initial  conditions, 
micrometeoroid  Impacts,  and  thermal  bending.  In  particular,  thermal 
bending  effects  are  so  serious  that  it  is  unlikely  that  the  present 
design  of  PACSAT  can  even  approach  the  straightness  requirements. 

A  partial  study  of  PACSAT  dynamics  has  been  conducted  by  Burke [11], 
who  considered  the  fundamental  mode  problem,  but  did  not  treat  the 
various  forces.  Some  of  Burke's  ideas  are  included  in  this  report. 


The  several  topics  are  investigated  in  the  following  sections.  The 
dynamic  problem  is  formulated  in  general  in  Section  11.  Several  mutually 
consistent  assumptions  are  employed  to  simplify  the  equations.  The 
fundamental  modes  of  the  simplest  problem'-circular  orbit  and  only 
gravity  gradient  and  elastic  forces--are  found  in  Section  III.  This 
proves  to  be  a  mathematical  problem  involving  singular  perturbations  and 
rather  delicate  analysis.  All  modes  have  real  frequencies,  so  the 
unperturbed  problem  is  stable  for  small  motions  (an  assumption  used  in 
the  analysis) .  The  effects  of  initial  conditions  are  evaluated  at  this 
point . 

Ellipticity  is  treated  in  Section  IV.  It  is  shown  that  small 
ellipticlty  drives  the  In-plane  motion  to  produce  a  libratlon  at  the 
orbital  frequency  with  amplitude  equal  to  the  eccentricity.  The 
frequencies  of  the  librational  modes  are  changed  slightly.  The  out- 
of -plane  equation  is  proved  to  be  unstable,  but  the  time  for  instability 
effects  to  appear  is  over  4,500  years.  There  is  no  significant  effect  on 
vibration. 

Earth  oblateness  effects  are  studied  in  Section  V.  Coupling  between 
equations  excites  a  librational  resonance  in  the  out-of-plane  motion  for 
near-circular,  near-equatorial  orbits,  but  the  amplitude  is  small,  being 
proportional  to  the  product  of  eccentricity  and  inclination.  In  Section 
VI,  solar  and  lunar  gravity  and  solar  radiation  pressure  are 
Investigated,  and  all  prove  to  have  small  effects. 

Micrometeoroid  impacts  are  studied  in  Section  VII.  Assuming  the 
particle  is  small  enough  that  it  does  not  destroy  the  integrity  of 
PACSAT,  the  effect  of  micrometeoroid  impact  is  closely  related  to  the 
effect  of  initial  conditions,  since  the  impact  is  equivalent  to  a  point 
discontinuity  in  velocity.  It  is  shown  that  the  effect  on  straightness 
may  be  quite  significant,  since  a  meteoric  particle  with  a  few  micrograms 
mass  can  produce  vibratory  motion  comparable  to  the  requirements,  and 
such  particles  may  be  encountered  as  often  as  once  a  month. 

Thermal  bending  effects  are  considered  In  Section  VIII.  The  solar 
heat  input  causes  differential  expansion,  which  makes  PACSAT  curve  away 
from  the  sun.  Eclipsing  by  the  earth,  or  shadowing  within  PACSAT,  causes 
the  solar  input  to  vary  in  time  to  excite  the  lowest  vibrational  mode  of 
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PACSAT  as  a  resonance,  and  the  resulting  displacement  is  major 
the  absorptivity  of  the  material  can  be  reduced  below  .0007,  a 
difficult  task,  the  straightness  requirements  cannot  be  met. 
Finally,  conclusions  appear  in  Section  IX, 


Unless 

very 
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II.  FORMULATION  OF  THE  GENERAL  PROBLEM 


The  geometrical  configuration  Is  shown  in  Fig.  1.  PACSAT  is  in  a 
near-circular,  near-equatorial  orbit.  It  is  nominally  vertical  and 
straight,  but  is  actually  somewhat  tilted  and  disturbed  into  a  sinuous 
shape.  The  figure  is  not  to  scale,  since  PACSAT  is  short  compared  to  the 
orbit  radius  and  the  distortions  are  small.  The  vector  r  is  the  radius 
vector  from  the  center  of  the  earth  to  a  point  P  on  PACSAT,  and  the 
directed  length  s  is  the  distance  from  the  center  of  mass  C  to  the  point 
on  the  undistorted  axis  of  PACSAT  corresponding  to  P.  With  L  the  length 
of  PACSAT  (1,500  m) ,  the  range  of  s  is  from  -L/2  to  L/2. 

A  series  of  assumptions  is  made  to  permit  analytic  treatment  of  the 
problem: 


Fig.  1  -  PACSAT  in  circular  equatorial  orbit 
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1.  Since  the  ratio  of  the  length  of  PACSAT  (1,500  m)  to  the  radius 

1  -4  I 

of  the  orbit  ("13,000  km)|  is  on  the  order  of  10  ,  all 

expressions  will  be  expanded  in  powers  of  this  ratio,  and  only- 

first-order  terms  kept. 

2.  All  slopes  and  curvatures  are  small,  so  in  forming  the 
Lagrangian,  only  second-order  terms  will  be  kept.  This  is 
consistent  with  the  requirement  of  small  displacements 
established  in  Section  I,  and  leads  to  linear  differential 
equat ions . 

3.  All  forces  other  than  gravity  and  elastic  forces  may  be  treated 
as  perturbations.  This  assumption  will  be  tested  as  each  force 
is  considered. 

4.  The  various  effects  will  be  considered  mutually  Independent, 
iinless  they  can  interact  to  produce  resonances. 

5.  The  PACSAT  structure  is  inextensible.  This  is  consistent  with 
the  properties  of  steel.  A  consequence  of  this  assumption  is 
that  there  are  no  longitudinal  vibrations . 

6.  The  PACSAT  rod  structure  is  very  long  compared  with  its 
diameter.  Hence,  the  axial  moments  of  inertia  may  be  neglected. 
Further,  there  will  be  no  torsional  vibrations. 

7.  The  PACSAT  structure  is  uniform,  and  the  motion  is  continuous. 

8.  There  is  no  damping  present.  Although  PACSAT  contains  "sticky 
hinges,"  the  gravity-gradient  tension  will  cause  them  to  lock. 

It  is  easily  shown  that  the  internal  friction  in  steel  is  not 
sufficient  to  provide  appreciable  damping. 

Under  these  assumptions,  the  equations  of  the  system  may  be 
formulated  and  simplified.  The  total  motion  may  be  separated  into  three 
parts,  the  motion  of  the  center  of  mass,  a  rigid-body  libration  about  the 
center  of  mass,  and  vibrations  transverse  with  respect  to  the  libration. 
The  formulation  will  use  the  Lagrangian  technique.  Although  other 
methods  could  have  been  employed,  the  Lagrangian  is  the  most  convenient 
when  thermal  bending  is  treated. 
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The  structure  of  PACSAT  in  the  undistorted  state  is  depicted  in 
Fig.  2.  The  central  steel  rod,  with  a  diameter  of  .132  cm,  carries  80,001 
spheres.  These  are  spaced  1.875  cm  between  centers  (1/2  wavelength  at 
3.75  cm),  and  have  radius  .595  cm  (3.75/2^),  the  radius  at  which  the 
radar  scattering  of  the  sphere  is  a  maximum.  The  spheres  may  be  solid 
aluminum,  through-drilled,  or  they  may  be  plastic,  to  save  weight,  coated 
with  aluminum  to  provide  radar  reflectivity. 

Following  Burke[ll],  we  assume  that  the  massive  solid  spheres 
prevent  bending  of  the  rod  within  the  sphere.  Hence,  the  effective 
length  of  the  rod  bearing  solid  spheres  is  taken  to  be  .685  x  1,500/1.875 
=  548  m.  The  plastic  spheres  are  sufficiently  light  that  they  do  not 
inhibit  the  bending,  and  the  effective  length  is  the  full  1,500  ra.  The 
effective  length  will  be  further  modified  when  thermal  bending  is 
considered.  This  use  of  effective  length  permits  us  to  treat  the  motion 
as  continuous,  as  per  assumption  7,  instead  of  inserting  immobile 
sections  and  having  to  investigate  a  very  complicated  boundary  value 
problem. 


Repeated  801X)1  times 


Fig.  2  —  PACSAT  structure 
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The  differential  equations  and  boundary  conditions  which  govern  the 
motion  will  be  derived  from  the  basic  Hamilton's  principle [ 12] , 


dt[T  -  V]  =  0 


2.1 


where  6  denotes  the  operation  of  variation,  T  is  the  kinetic  energy,  and 
V  is  the  potential  and  strain  energy.  The  energies  must  be  expressed  in 
terms  of  the  coordinates.  As  shown  in  Fig.  1,  the  position  of  a  point  P 
is  described  by  a  vector  r  which  determines  the  distance  of  P  from  the 
center  of  the  earth,  and  by  a  distance  s  which  identifies  the  point  P. 
Under  the  assumption  that  the  rod  is  Long  compared  with  its  radius,  the 
internal  distortions  of  the  cross  section  of  the  rod  will  be  neglected, 
and  it  will  be  assumed  that  the  position  of  each  point  on  the  cross 
section  is  characterized  by  the  distance  s.  This  also  corresponds  to  the 
assumption  of  small  displacements,  slopes,  and  curvatures.  Under  these 
assumptions,  the  motion  of  the  rod  may  be  described  by  writing  the 
position  vector  r  as  a  function  of  s  and  t,  the  time,  and  ignoring  the 
internal  coordinates. 

Let  p  denote  the  mass  density  per  unit  length  of  the  structure.  For 
the  steel  rod  with  alxuninum  spheres,  the  total  mass  per  1.875  cm  section 
is  the  mass  of  the  rod  part  plus  the  mass  of  the  sphere.  The  volume 
density  of  steel  is  7.86  g/cm^,  that  of  aluminum  is  2.69.^  The  mass  of 
the  steel  rod,  radius  .066  cm,  length  1.875  cm,  is  .202  g.  The  volume  of 
the  entire  sphere  (radius  .595  cm)  is  .882  cm  ,  that  of  the  section  of 
rod  within  the  sphere  is  .0163  cm,  and  hence  the  mass  of  the  drilled 
aluminum  sphere  is  2.330  g,  the  total  mass  per  1,875  cm  section  is 
2.531  g,  and  the  linear  mass  density  p  is  1.35  g/cm.  For  the  plastic 
spheres  loaded  rod,  we  assume  that  the  spheres  are  sufficiently  light  that 

^ These  and  all  other  numerical  values  pertaining  to  materials  are 
taken  from  International  Critical  Tables ,  McGraw  Hill  Book  Co. ,  New  York, 
1927. 
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the  mass  density  is  that  of  the  rod  alone,  which  from  the  above  figures 
is  .107  g/cm.  These  values  were  first  derived  by  Burk.e[ll]. 

From  here  on,  PACSAT  will  be  referred  to  as  "the  rod,"  except  where 
the  spheres  are  to  be  specifically  discussed.  The  total  kinetic  energy 
of  the  rod  is 


2 


All  integrations  over  time  in  the  ensuing  derivation  have  the  limits  t^^, 
t^;  all  integrations  over  s  have  the  limits  -L/2,  L/2.  The  limits  will  be 
omitted  and  implicitly  understood.  According  to  the  principle  of 
variation  as  discussed  in  Ref.  12,  the  time  integral  of  the  Lagrangian  is 
stationary  with  respect  to  variations  of  the  coordinates  about  the 
"trajectory"  which  connects  the  state  of  the  system  at  the  initial  time 
tj^  to  the  state  at  the  final  time  t^,  provided  these  variations  vanish  at 
the  endpoints.  Then  the  variation  of  the  kinetic  act^n  is 


6T 


=  p/dt/ds-^ 


2.3 


The  definition  of  variations  permits  the  commutation  of  the 
operations  of  variation  and  time  differentiation.  Integrating  by  parts, 


6r 


2-> 

^  -  /dt— 
3t 


6T  =  p/ds 


2.4 
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The  integrated  term  vanishes,  since  the  variations  must  vanish  at  the 
endpoints.  The  variation  of  the  kinetic  action  is  thus 


5T  =  -p/dt/ds- 


^2-»- 
d  r 


3t 


6r 


2.5 


and  is  immediately  recognizable  as  the  acceleration  term. 

The  potential  energy  is  composed  of  the  gravitational  potential  and 
the  work  done  by  the  external  forces.  In  most  cases  the  external  forces 
are  derivable  from  a  potential,  but  the  more  general  case  can  be  treated 
just  as  easily.  Let  U(r)  be  the  gravitational  potential  per  unit  length: 


U(r)  =  -pv/|rl 


where  ji  is  the  product  of  Newton's  constant  and  the  mass  of  the  earth, 
and  let  F(r,t)  denote  the  external  forces  per  unit  length.  The  potential 
energy  associated  with  these  external  forces  is 


W(r,t)  =  -jT  •  (r'  -  rg)dr’ 


2.7 


where  r^  is  an  initial  state  of  the  system.  The  variation  of  the 
gravitational  action  is 


6U  =  -6fdt/dspii/|r(s,t)  I  =  t/dtldspuCr  •  6r)/|r| 


2.8 
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and  the  variation  of  the  external  potential  action  is 


6W  =  -/dt/dsF(r,t)  •  6r  . 


2.9 


The  strain  energy  is  found  by  the  methods  of  the  theory  of 
elasticity.  The  strain  in  a  rod  element  may  be  resolved  into  a  tension 
along  the  rod  and  a  bending  couple  around  the  rod  axis.  For  the  tension, 
see  Fig.  3a.  The  action  of  the  tension  is  to  stretch  the  rod  element, 
originally  of  length  ds ,  into  a  section  of  length  d£.  The  work  done  is 
the  force  times  the  elongation.  The  change  in  rod  length  is 
approximately  l/2(3r/3s)^ds , ^  and  the  strain  energy  of  the  entire  rod  due 
to  tension  is 


Fig.  3a  —  Tension  in  a  rod  element 


^Here  and  hereafter  the  square  of  a  vector  quantity  means  the  square 
of  the  magnitude  of  the  vector. 
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V 


T 


j/dsT(s,t) 


2,10 


Varying  the  time  integral  of  this  quantity, 


6/dtV^  =  /dt/dsT(s,t)|j 


2.11 


The  variation  operation  and  length  differentiation  can  be  commuted. 
Integrating  by  parts. 


6/dtV^ 


2.12 


Since  the  variations  need  not  be  constrained  at  the  ends,  this 
imposes  the  boundary  condition: 


T(s,t)  =0  s  =  ±L/2  . 


2.13 


This  condition,  that  the  rod  cannot  sustain  tension  at  its  ends,  is 
natural  to  the  physics  of  the  problem. 

For  the  bending,  see  Fig.  3b.  Suppose  that  the  initially 
undistorted  rod  is  bent  corresponding  to  a  radius  of  curvature 
Consider  the  extension  of  a  filament  of  the  rod  at  a  distance  y  from  the 
center  of  the  section.  The  length  of  the  filament  is  altered  by  the 
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Fig.  3b  -  Bending  of  a  rod  element 


L 


bending  in  the  ratio  1/(1  +  y/R^) .  Thus,  on  the  side  of  the  axis 

for  which  y  is  positive,  the  outward  side,  the  filament  is  extended, 

while  on  the  inward  side  it  is  compressed.  By  the  definition[ 13]  of 

Young's  modulus,  E,  the  force  necessary  to  produce  the  extension  y/R  over 

the  element  of  area  dA  is  EdAy/R„,  The  couple  by  which  the  bending  is 

K 

resisted  amounts  to 


E/dAy  '  y/R^  =  Zl/R^ 


2.14 


where  I  is  the  moment  of  inertia  of  the  area  (not  the  mass)  around  the 
bending  axis.  The  angle  of  bending  corresponding  to  a  length  of  axis  ds 
is  ds/R^.  Since  the  mean  value  is  half  the  final  value  of  the  couple, 
the  work  required  to  bend  ds  to  curvature  is 
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~EIds/R^  . 


2.15 


This  argument  is  substantially  due  to  Lord  Rayleigh[ 14] . 

The  radius  of  curvature  is  the  reciprocal  of  the  curvature  itself. 
By  the  assumption  of  small  slopes,  the  curvature  can  be  replaced  by 
3  r/3s  .  Hence,  the  strain  energy  of  the  rod  due  to  bending  is 


’'b  - 


3s' 


2.16 


Forming  the  variation  of  the  time  integral  of  Vg, 


6/dtV„ 


=  El/dt/ds 


2“^ 
.3  r 


3s 


2.17 


Again  the  variation  and  s  derivative  can  be  commuted.  Integrating  by 
parts  twice. 


6/dtV„  =  El/dt 

Q 


n 

33^ 


3s 


6r  - 


3s^ 


6r 


+  /ds 


.L 

2 


3s^ 


St 


2.18 


Since  the  variation  6r  and  its  s  derivative  are  independent  arbitrary 
quantities,  this  yields  two  further  boundary  conditions.  Finally, 
collect  all  the  terms  proportional  to  6r  under  the  Integral  sign  in  the 
expressions  2.5,  2.8,  2.9,  2.12,  and  2.18.  The  arbitrariness  of  6r  along 
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the  rod  requires  that  its  coefficient  vanish,  yielding  the  governing 
differential  equation.  Thus,  the  problem  to  be  solved  is  formulated  as 
the  differential  equation: 


9^? 

^  2 
at 

-pnr 

3 

r 

4^ 

-Elif. 

as 

F(r,t)  . 

2.19 

inertia 

Gravity 

Twwfen 

Bent^ng 

Exterrnri 

The  boundary  conditions  are 


At  s  =  tL/2  T(s,t)  =  0 


2.20a 


2-^  2 
a  r/3s 


0 


2.20b 


a^r/Ss^ 


0 


2.20c 


These  mean  that  at  the  ends  the  rod  cannot  sustain  tension,  bending,  or 
shearing,  all  of  which  are  clearly  required  conditions. 

In  deriving  2.19  and  2.20,  we  have  assumed  that  the  kinetic  energy 
is  unaffected  by  the  bending.  Actually,  there  is  kinetic  energy 
associated  with  the  rotatory  motion  of  the  rod,  the  so-called  "rotary 
inertia"  term(14j,  but  this  contributes  a  term  to  the  differential 
equation  which  is  of  the  order  of  magnitude  of  the  square  of  the  ratio  of 
the  rod  radius  to  its  length,  relative  to  the  potential  energy  term 
derived  from  2.18,  and  this  can  safely  be  ignored.  This  neglect  of  the 
"rotary  inertia"  is  customary  in  vibration  problems. 

The  condition  of  in extensibility  leads  to  a  restriction  on  the 
possible  motions.  Consider  two  points  on  the  rod,  which  in  undistorted 
motion  are  at  positions  s  ,  s  ,  separated  by  a  distance  ds .  As  the  rod 
bends,  they  will  be  at  positions  r(s^,t),  rCs^jt).  The  distance  between 
the  points  must  remain  unchanged,  so, 
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|r(s^  +  ds,t)  -  r(s^,t)|  =  ds  . 


2.21 


Dividing  by  ds  and  taking  the  limit  as  ds  approaches  zero  yields 


9r(s,t)  ^  ^  2  22 

3s  * 


This  relation  will  be  used  for  simplification  of  the  equations. 

We  first  show  that  equations  2.19  and  2,20  imply  the  conservation  of 
angular  momentum.  This  is  in  contradistinction  to  the  analysis  of 
Burkefll],  who  imposes  angular  momentum  conservation  as  an  additional 
constraint  and  deduces  that  the  librational  motion  cannot  exist,  a  result 
clearly  in  disagreement  with  observation.  To  show  conservation  of 
angular  momentum,  take  the  vector  cross  product  of  the  vector  r  and  the 
vector  differential  equation  2.19,  then  Integrate  over  the  length  of  the 
rod.  For  the  left  side,  we  have  (limits  of  integration  -L/2,  L/2 
implied) 


p/dsr  X 


3 


2-^ 

r 


d  rj 

^p/dsr 


2.23 


This  expression  is  immediately  recognized  as  the  time  derivative  of  the 
total  angular  momentum  of  the  rod.  The  first  term  on  the  right  is  zero 
(r  X  r  =  0) .  For  the  second  term,  integrate  by  parts 


fdsT  X  |_T(s,t)|| 


=  r 


L 

2 

_L 

2 


;dsT(»,t)|f  .  2.24 
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The  first  term  on  the  right  side  of  2.24  vanishes  by  the  boundary 
condition  2.20a,  the  second  from  the  cross  product  of  identical  vectors. 
For  the  third  term  on  the  right,  integrate  by  parts  twice: 


fdsr  X 


d  r 

a  ^ 
9s 


ds^ 


3s 


3s^/ 

2 


9  2? 

+  ids — j  X 
9s^ 


9s^ 


2,25 


The  integrated  terms  vanish  by  the  boundary  conditions  2.20b  and  2.20c, 
the  Integral  from  the  cross  product  of  identical  vectors.  Hence,  the 
result  of  integration  is 


■^pfdsr  ^  ^  =  /dsr  x  F(r,t) 


2.26 


Equation  2.26  asserts  that  the  rate  of  change  of  angular  momentum  equals 
the  torque  applied  by  the  external  forces.  If  there  are  no  external 
forces,  the  angular  momentum  is  conserved.  Hence,  no  additional 
constraints  need  be  applied  to  2.19  and  2.20,  and  Burke's  result [11}  is 
erroneous . 

The  differential  equation  2.19  and  boundary  conditions  2.20  describe 
the  motion  of  a  point  on  the  rod,  relative  to  the  center  of  the  earth. 
However,  PACSAT  is  known  to  be  in  orbit,  executing  small  motions  around 
the  local  vertical.  Hence,  we  separate  the  motion  into  the  orbital, 
librational,  and  vibrational  motions,  and  treat  the  several  motions 
separately.  To  this  purpose,  write 

r(s,t)  =  R(t)  +  su(t)  +  w(s,t)  =  R(t)  +  v(s,t)  .  2.27 

Center  Center 

of  mass  Libration  Vibration  of  mass  Displacement 
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The  vector  R(t)  represents  the  motion  of  the  center  of  mass  of 
PACSAT.  It  is  of  the  order  of  magnitude  of  the  distance  of  PACSAT  from 
the  center  of  the  earth  (1.2  -  1.5  x  10  km).  The  expression  su(t)  is 
the  rigid-body  libration  about  the  vertical  (the  displacement  is 
proportional  to  the  distance  from  the  center  of  mass),  and  is  of  the 
order  of  the  length  of  PACSAT  (1.5  km).  The  vector  w(s,t)  measures  the 
distortion  of  PACSAT  from  straightness.  This  distortion  executes  free 
and  forced  vibrations  about  the  librating  rigid  body.  It  is  of  the  order 
of  magnitude  of  the  permissible  displacement  (.5  cm).  Hence,  the  terms 
of  2.27  are  such  that  expansion  of  the  differential  equation  2.19, 
boundary  conditions  2.20,  and  subsidiary  condition  2.22  in  ratios  of  the 
libration  to  center  of  mass,  or  vibration  to  libration,  are  all  rapidly 
convergent,  and  only  first-order  terms  need  be  kept.  We  further  assume 
that  the  scale  length  of  variation  of  the  perturbing  force  F  is  at  least 
on  the  order  of  magnitude  of  R.  This  is  satisfied  for  the  oblateness 
(order  R) ,  lunar  gravity  (order  distance  earth  to  moon  =4x10  km) ,  and 
solar  gravity  and  radiation  pressure  (order  distance  earth  to  sun 

g 

=  1,5  X  10  km),  all  the  forces  we  will  consider. 

The  tension  T(s,t)  may  also  be  expanded.  As  is  well  known[9],  a 
rigid,  vertical,  straight  rod  in  circular  orbit  experiences  a  tension 
which  is  parabolic  along  its  length.  Hence,  we  take  for  the  tension  the 
form: 


T(s,t) 


=  T^Ct) 


+  T^(s,t) 


2.28 


where  the  first  term  on  the  right  represents  the  tension  associated  with 
the  center-of -mass  and  libratlonal  motion,  displaying  a  parabolic 
variation  which  vanishes  at  the  ends  (2.20a);  the  second  term  is  the 
additional  tension  associated  with  vibration. 
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The  gravity  term  expands  as  follows: 


r 

3 

r 


_ R  +  V 

+  2R  •  V  +  V 


2^3/2 


(R  +  v) 


2.29a 


2.29b 


*  8 

where  the  neglected  terms  are  on  the  order  of  10  relative  to  the 
leading  term.  The  external  force  expands  in  the  form: 


FCr,t)  s  T(R,t)  +  [VF(R,t)]  •  v 


2.30 


where  the  expression  VF(R,t)  denotes  the  3x3  dyadic  which  in 
rectangular  coordinates  has  the  representation: 


VF(R,t) 


aF 

8F 

3F 

X 

X 

X 

8X 

9Y 

3Z 

3F 

3F 

3F 

y 

—L 

3X 

3Y 

3Z 

3F 

3F 

3F 

z 

z 

z 

ax 

3Y 

3Z 

2.31 
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and  which  transforms  under  coordinate  transformations  like  a  tensor  of 
second  rank.  All  components  of'2.31  are  evaluated  at  the  point  R.  All 
other  terras  in  2.19  and  2.20  are  linear  and  have  obvious  expansions.  We 
obtain  for  the  center  of  mass  motion: 


+  ^F{R,t) 
R  ^ 


2.32 


for  the  libration; 


8T„(t)^  1 

— — — u  +  ^(VF) 
t2  P 

pL 


u  , 


2.33 


and  for  the  vibrational  motion: 


-h 

w 


u 


2.34 


s  34  +  i(,|) 

»  a,*  '■ 


w 


The  boundary  conditions  become 


"2 


=  0 


At  s  = 


2.35 
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The  expansion  of  the  condition  2.2  yields 


|u(t)|^  =  1 


u  •  w  =  0 


jw|  negligible  . 


2.36a 


2.36b 


2,36c 


The  condition  2.36a,  that  the  vector  u  have  unit  magnitude,  demonstrates 
the  rigid-body  character  of  the  libration.  The  condition  2.36b  shows 
that  the  vibrational  motion  is  always  transverse  to  the  libration,  and 
2.36c  demonstrates  that  higher  order  terras  in  the  vibration  are  to  be 
neglected.  Thus,  we  have  the  motion  of  the  center  of  mass  under  central 
gravity  and  external  forces,  as  described  by  2.32,  the  rigid  libration 
about  the  center  of  mass  under  gravity  gradient,  tension,  and  external 
forces,  described  by  2.33  and  2.36,  and  the  linearized  transverse 
vibration  about  the  libration  under  gravity  gradient,  tension,  elastic 
bending,  and  external  forces,  with  no  bending  or  shearing  at  the  ends, 
described  by  2.34,  2.35,  and  2.36b.  These  equations,  under  the 
appropriate  restrictions,  will  be  considered  in  the  subsequent  sections. 
The  equations  will  be  modified  when  thermal  bending  is  considered. 

The  assumption  of  linearity  may  be  relaxed  if  we  are  willing  to 
accept  a  complicated  computer  program  and  implement  numerical  methods  for 
treating  the  nonlinear  equations.  A  recent  analysis  by  Kane  and 
LevinsonllS,  16]  describes  such  a  treatment,  but  the  limitations  on  the 
motion  required  to  meet  the  requirements  on  radar  cross  section  make  such 
complications  unnecessary. 
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III.  THE  BASIC  PROBLEM 


In  this  section,  the  simplest  problem  will  be  considered.  The  orbit 
of  PACSAT  will  be  taken  as  circular  and  equatorial.  All  perturbing 
forces  will  be  neglected.  We  shall  find  the  modes  of  libration  and 
vibration  of  PACSAT  under  these  circumstances,  and  determine  the  effects 
of  initial  conditions. 

The  motion  is  conveniently  described  in  cylindrical  coordinates. 
Referring  to  Fig.  1,  the  distance  R  to  the  center  of  mass  of  the 
satellite  will  be  constant.  The  angle  0,  which  locates  the  position  of 
the  center  of  mass  in  the  orbital  plane,  measured  from  a  reference 
direction  (the  first  point  of  Aries  is  convenient),  will  Increase 

9 

linearly  with  time  at  the  orbital  rate  0,  which  is  also  constant.  The 
displacements  u  and  w  will  be  expressed  in  terms  of  i^,  ig ,  and  i^,  the 
unit  vectors  in  the  radial,  velocity,  and  normal  to  the  orbit  plane 
directions.  If  we  also  use  rectangular  coordinates,  with  the  z  axis  in 
the  direction  of  the  North  Pole  and  the  x  direction  toward  Aries,  we  have 
the  time  dependences 

r,,  =  1  cosSt  +  1  sin0t  3.1a 

R  X  y 

1-  =  -i  sin0t  +  i  cos6t  .  3.1b 

0  X  y 


->■  -*■ 
The  vector  R  of  equation  2.32  is  simply  Rig- 

are  expressed  as 


The  vectors  u  and  w 


^  ^  ^2^0  “3^2 

w  =  w +  w^ig  +  W3i^  . 


3.2a 


3.2b 
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The  conditions  2.36  yield 


2  2  ^  2  T 

^1  +  '^2  ^3  ^ 


3.3a 


UiW^  +  U2W2  +  i^3W3  =  0  . 


3.3b 


If  the  rod  were  vertical  and  straight,  we  would  have  =  1,  all 
other  terms  zero.  The  condition  3.3a  then  permits  the  representation 


u 


1 


=  COSXCOS^ 


3 .4a 


U2  =  cosxsin^ 


3.4b 


=  sinX  . 


3.4c 


The  angle  ^  is  the  in-plane  angular  displacement  from  vertical,  the  angle 
X  is  the  out-of“plane  angular  displacement  from  vertical.  These  will  be 
treated  as  first-order  quantities.  The  displacements  and  w^  are  the 
transverse  displacements  from  straightness,  also  to  be  treated  as  first 
order.  Equation  3.3b  shows  that  the  radial  displacement  w^^  is  a  second- 
order  quantity. 

The  angular  velocity  0  and  orbit  radius  R  satisfy  the  relation 


R^0^  =  P  =  3.986  X  10^  km^/sec^  .  3.5 

A  corresponding  relation  C4.1e)  holds  in  elliptic  orbits.  For  PACSAT  at 
an  altitude  of  1  earth  radius  (6,378  km),  the  orbital  period  is  3  hours, 
59  minutes;  at  1.5  earth  radius,  the  period  is  5  hours,  34  minutes. 
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The  velocities  and  accelerations  must  be  represented  in  the 
cylindrical  coordinates.  Differentiating  3.1  yields 


Ir,  = 


lo  = 


It- 


irt 


•-i- 

3.6a 

-ei_ 

3.6b 

R 

•2"^ 

-e 

3.6c 

K 

®  0  ■ 

3.6d 

The  acceleration  u  is  given  by 


u  =  u^ij^  +  +  u^ij^  +  ... 


3.7a 


=  (u^  -  2eu2  -  **■  ^^2 


Ths  gravity ■“gra.dien’t  term  on  tlie  right  sidfe  of  2,33  becomes 


-0  (-2u^i^  +  u^ig  +  U3i^) 


3.8 


There  are  corresponding  terms  for  w. 

We  first  consider  the  libration.  With  the  perturbing  force 
omitted  and  some  simplification,  2.33  becomes 


.2  .. 

-  30  -  20u2  - - 2^1 

pL 


3.9a 
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20U, 


+  0 


8T| 

pL 

8T 

I 

pL 


0 

2'‘2 


0 

2’"3 


3.9b 


3.9c 


Multiply  the  first  equation  by  u^^,  the  second  by  U2,  the  third  by  u^,  and 
add.  This  provides  an  equation  for  the  tension  T^.  Equations  3.4  may  be 
differentiated  to  express  the  derivatives  of  u^,  u^,  and  in  terms  of 
Ip,  X,  and  their  derivatives.  It  is  convenient  to  introduce  the  angle  6 
as  the  independent  variable.  Derivatives  with  respect  to  B  will  be 
denoted  by  a  prime.  Also,  multiply  3,9a  by  u^,  3.9b  by  u^,  and  subtract. 
This  will  yield  an  equation  which  is  appropriate  for  determining  ip . 
Einally,  substitute  the  expression  for  the  tension  into  the  right  side  of 
3.9c  and  simplify.  This  equation  is  suitable  for  X.  After  the 
manipulations,  there  results: 


Tn  =  +  cos^x(l  +  ^’)^  +  3cos^Xcos^tp  -  1]  3.10a 

U  o 

+  Scos^sintp  =  2tanX(l  +  ip')x'  3.10b 

x"  +  [(1  +  'p')^  +  3cos^ip]s±nxcosx  =  0  .  3.10c 


To  this  point,  there  have  been  no  simplifications  produced  by  the 
small  values  of  ip  and  X.  Now  we  wish  to  restrict  the  libratlon  amplitude 
to  1.5  deg  =  .025  rad.  Keeping  terms  of  orders  0  and  1  in  equation  3.10a 
and  of  orders  1  and  2  in  equation  3.10b  and  c  should  be  more  than 
sufficient  accuracy.  Thus,  we  expand 


,p  =  ,p^  +  ip^ 


3.11a 
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X 


3.11b 


and  obtain  the  set  of  equations: 


Tq 

3.12a 

+  3^^-  0 

3.12b 

x'l  +  0 

3.12c 

+  3#2  =  2Xj^x| 

3.12d 

X2  +  ^X2  =  . 

3.12e 

The  first  equation,  3.12a,  represents  the  gravity-gradient  tension, 
which  has  been  derived  many  times  previously[9 ,  10,  and  other  papers],  as 
modified  by  libration.  For  PACSAT  at  1  earth  radius,  p  =  1.35  g/cro, 

L  =  1,500  m,  0  =  2ir/14,340  rad/sec,  the  tension  is 


Tp  =  2,187(1  +  dynes  . 


The  full  length  is  used,  since  the  entire  body  partakes  of  the  rigid 
motion.  This  value  was  also  obtained  by  Burke[ll]. 
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The  equations  3.12b  and  3.12c  solve  Immediately  as 


tpT  =  H — psin/To 

10  /3 


Xi  =  Xocos20  +  |xosiii20 


3.l4a 


3.14b 


the  subscript  0  denotes  initial  value  and  0  is  measured  from  the 
initial  time.  Thus,  the  first-order  in-plane  libration  oscillates  at  a 
frequency  equal  to  ^3  times  the  orbital  frequency,  the  first-order  out- 
of -plane  libration  oscillates  at  twice  the  orbital  frequency.  These 
results  also  are  well  known[9,  10].  If  we  assume  the  initial  libration 
velocities  are  zero,  which  should  be  at  least  approximately  satisfied  for 
any  reasonable  deployment  scheme,  the  sin  terras  may  be  dropped. 

The  variation  of  the  tension,  as  shown  by  equations  3.12a  or  3.13, 
has  a  maximum  equal  to  Z’pn/ if  the  initial  displacement  is  entirely  in 

U  Q 

the  in-plane  direction.  For  an  initial  value  1.5  ,  the  second  term  lu 
3.12a  or  3.13  is  at  most  2.9  percent,  so  the  tension  is  substantially 
constant . 

The  equations  3.12d  and  3.12e  may  be  solved  easily,  and  yield: 


3.15a 


X2  =  1  +  /3)0  +  (4  +  /3)sin(2  -  V^)  0  -  5sin20] 

3.15b 


where  the  second-order  terms  have  been  fitted  to  zero  initial  conditions. 
If  the  motion  is  limited  to  1.5  deg,  the  maximum  values  of  these 
expressions  are  respectively  .00032  (Xq  =  .0262),  and  .00033  (;(ip  — 

=  .0185).  Thus,  the  second-order  corrections  to  the  libration  are  on  the 
order  of  1  percent  and  may  be  neglected. 
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The  first-order  libration,  3.14,  represents  a  complicated  time 
dependence  for  the  orientation  of  the  libration.  The  motions  in  the  two 
planes  are  at  different  frequencies,  so  the  attitude  orientation  of  the 
body  follows  a  Lissajous  figure.  If  we  assume  the  initial  displacement 
is  of  unit  amplitude  and  displaced  45  deg  out  of  plane,  the  resulting 
motion  for  the  first  four  orbits  is  as  plotted  in  Fig,  4,  The  curve  may 
be  thought  of  as  the  motion  of  the  tip,  plotted  in  the  local  east  (in¬ 
plane)  and  north  (out-of -plane)  coordinate  system  (a  plane  perpendicular 


ln-pl«ne  relative  dispiscemsnt 


'  Fig.  4  —  Libration  Lissajous  pattern 

I 
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to  the  radius  vector) .  The  numbers  at  the  top  and  bottom  of  the  curve 

represent  the  time  in  fractions  of  an  orbit.  The  motion  begins  at  t  =  0 

in  the  upper  right  corner  of  Fig.  4,  and  is  at  the  top  each  half  and 

whole  orbit,  at  the  bottom  each  quarter  orbit,  as  can  be  seen  from  3.14b. 

Since  the  frequencies  in  the  two  planes  are  incommensurable,  the  curve 

will  never  repeat,  and  will  eventually  cover  the  entire  square.  However, 

it  can  come  very  close  to  repeating.  Table  2  shows  some  near  repetitions 

of  the  initial  value.  The  first  column  is  the  orbit  number,  the  second 

the  elapsed  time,  and  the  third  the  leftward  displacement  from  the  upper 

2  2 

right  corner  in  imits  of  the  initial  displacement  +  X^)  =  1. 

The  closeness  to  repetition  at  such  relatively  low  values  of  orbit 
number  is  remarkable.  The  34.69  day  value  corresponds  to  the 
approximation 


/3  =  1.7320508  ~  362/209  =  1.7320574  . 


3.16 


Closeness  to  other  corners  of  the  square  may  be  even  better,  but  does  not 
represent  pattern  repetition. 


Table  2 


LISSAJOUS  REPETITIONS 


Orbit  Number 

Elapsed  Time  (days) 

Displacement 

4 

.664 

.070747 

11 

1.83 

.038208 

15 

2.49 

.005159 

41 

6.80 

.002766 

56 

9.29 

.000371 

153 

25.39 

.000199 

209 

34.69 

.000027 

571 

94.77 

.000014 

780 

129.46 

.000002 
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It  is  shown  in  Ref.  3  that  the  mutual  viewing  time  o£  PACSAT  in 
equatorial  orbit  of  an  altitude  of  one  earth  radius  is  about  45  minutes 
for  two  observers  at  the  NE  and  NW  corners  of  the  continental  United 
States,  the  minimum  time,  and  is  about  110  minutes  for  adjacent  observers 
at  SO^N,  the  maximum  time.  These  correspond  to  .188  and  .460  part  of  an 
orbit.  In  the  first  case,  about  three-quarters  of  an  arc  of  Fig.  4  will 
be  covered,  in  the  second  case  nearly  two  arcs  will  be  traversed.  Thus, 
practically  the  full  range  of  libration  orientations  will  be  presented  to 
the  observers,  and  the  attendant  reductions  in  cross  section  will  have  to 
be  overcome . 

This  completes  our  discussion  of  free  libration,  driven  only  by 
initial  conditions.  Before  discussing  the  vibrational  motion,  we  briefly 
digress  to  consider  the  accuracy  of  the  assumption  that  the  rod  is 
inextensible.  With  the  possible  exception  of  forces  induced  by 
micrometeoroid  impacts,  which  we  will  treat  at  the  appropriate  occasion, 
the  gravity-gradient  tension  is  the  strongest  force  acting  to  produce 
extension  of  the  rod. 

The  total  elongation  of  the  rod,  obtained  by  integrating  the  tension 
over  the  length,  then  dividing  by  the  product  of  the  cross-sectional  area 
and  by  Young's  modulus  (2  x  10  dynes/cm  for  steel),  is  given  by 


As 


0  T  L 
_  ° 

3  ^ 


3.17 


For  PACSAT,  this  has  the  value  .03  cm.  Considering  the  total 
manufacturing  tolerances  and  possible  deformations  during  deployment, 
this  is  certainly  not  a  significant  quantity.  The  stretching  of  a 
section  of  length  1.875  cm  (sphere -to- sphere  spacing)  at  the  center  is 
6  X  lo”^  cm,  or  60  Angstrom  units,  far  below  any  manufacturing  tolerance. 
The  assumption  of  inextensibility  is  certainly  verified. 

We  DOW  take  up  the  subject  of  free  vibratory  modes.  This  involves 
solving  equation  2.34  with  the  boundary  condition  2.35  and  subsidiary 
condition  2.36b.  Only  first-order  quantities  will  be  kept,  which  implies 
that  the  effect  of  the  libration  on  the  vibratory  motion  will  be 
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neglected.  We  have  carried  out  a  more  detailed  analysis,  and  shown  that 
there  are  no  first-order  modifications  to  the  frequencies,  but  there  may 
be  second-order  effects.  Since  the  libration  is  limited  to  .025,  these 
are  changes  of  less  than  0.1  percent  and  will  not  be  considered  further. 

the  representation  3.2b  for  w  is  substituted  into  equation 
2.34,  the  only  surviving  first-order  terms  in  the  radial  component 
determine  the  correction  to  the  tension  as 


,2 

—  =  -2pe  w- 


3.18 


Anticipating  later  results,  we  have  the  order-of -magnitude  relation 


~  w/L  =  5 


10'^/1,500 


3.3  X  10 


3.19 


so  the  change  in  tension  is  very  small. 

The  in-plane  and  out-of-plane  components  of  the  vibratory  motion 
yield  the  differential  equations 


with  boundary  conditions 
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At 


2  3 

3  W2  3 

^  ~  3~ 

3s  3s 


0 


^2  .3 

3  3 

i~  ^  3~ 

3s  3s 


3.21a 


3.21b 


We  shall  use  8  as  one  independent  variable,  and  introduce  the  normalized 
length  X  as  the  other  where 


s 


-1  <  X  <  1  . 


3.22 


In  addition,  we  define  the  elasticity  parameter  B  as 


B  =  32EI/3pL^0^  . 


3.23 


The  parameter  B  measures  the  relative  effects  of  the  elasticity  and 
gravity-gradient  tension  on  the  motion.  For  PACSAT,  at  one  earth  radius 
altitude,  we  have  the  values 


1  ?  2 

E  =  2  X  10  dynes/cm  3.24a 

I  =  |■(,066)'^  =  1.49  X  lo"^  cm^  3.24b 


0  =  2Tr/14,340  sec 


3.24c 


For  the  two  cases  of  bare  rod  and  rod  loaded  with  aluminum  spheres,  we 
have  respectively  p  =  .107  g/cro,  L  =  1,500  m,  and  p  =  1.35  g/cm, 

L  =  548  m.  These  yield 


e  =  1.36  X  10 


rod  with  spheres  . 


3.25b 


4 


Thus,  p  is  a  small  quantity. 

We  see  that  the  differential  equations  and  boundary  conditions  are 
uncoupled,  so  in-plane  and  out-of-plane  motions  vibrate  Independently. 
Suppose  the  characteristic  frequencies  of  the  in-plane  (out-of -plane) 
motion  are  denoted  by  Mjp  .  Then  it  is  apparent  from  the  left 

sides  of  3.20  that 


2 

‘^0P,k 


IP,k 


+  0 


3.26 


Hence,  only  the  in-plane  motion  need  be  treated  in  detail.  The  out- 
of-plane  motion  will  be  deduced  by  comparison  of  coefficients  and  the  use 
of  3.26.  We  note  that  the  libration  frequencies  (Wjp^O  “  “OP,0  ~ 

also  satisfy  3.26,  from  which  we  are  justified  in  calling  the  libration 
the  lowest  odd  (in  z)  vibration  mode. 

With  the  new  independent  variables,  we  have 

3.27a 

3.27b 

These  may  be  solved  by  separation  of  variables,  thus 


2  „  _  3 


(  2. ^”2 


4 

3  w„ 


-  &— 


9x 


At  X  =  ±1 


2  3 

9  ^  ^2 


=  0  . 


9x 


9x' 


oo 

W2  =  SAj^(9)X^(x) 


3.28 
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which  in  turn  yield  the  equations: 


's;  +  "k 


2  2.  = 
k\  dx 


'3“ 


n 

(1  -  X 


2  3 

d  K  d  X 

At  X  =  ±1  ^  =  0  , 

dx  dx 


3.29 


3.30a 


3.30h 


The  solution  of  3.29  is 


A^(0) 

A,  =  A,  (O)cosu,  0  +  - sino),  9  . 

X  k  k  to,  k 

k 


3.31 


The  principal  task  of  this  section  is  to  investigate  the  eigenvalue 
problems  described  by  3.30. 

With  some  changes  of  notation,  this  is  the  same  as  the  eigenvalue 
problem  studied  by  Burke[ll].  He  actually  investigated  a  set  of  problems 
which  approximate  3.30,  but  his  treatment  of  the  full  problem  did  not 
consider  the  fact  that,  although  the  elastic  term  has  a  small 
coefficient,  the  problem  is  of  the  character  of  a  singular  perturbation, 
and  limit  procedures  cannot  be  applied.  Therefore,  Burke  s  results, 
which  are  very  limited  because  he  considered  only  the  lowest  eigenvalue, 
must  be  regarded  at  best  as  suspect. 

The  singular  nature  of  the  problem  arises  from  the  presence  of  the 
small  coefficient  in  the  highest  derivative  term  in  3.30a.  Suppose  P 
approaches  zero.  For  P  =  0,  the  differential  equation  3.30a  reduces  to 
Legendre's  equation.  The  eigenvalues  are  given  by 


|k(k  +  1) 


= 


1/2 
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and  the  eigenfunctions  are  the  Legendre  polynomials  Pj^(x).  But  this 
solution  does  not  satisfy  the  boundary  conditions  3,30b.  Hence,  the 
limit  of  the  solution  of  3.30  as  B  approaches  zero  is  not  the  solution 
for  B  =  0,  displaying  the  singular  character.  We  expect  that  any  attempt 
to  apply  a  perturbation  theory  corresponding  to  an  expansion  in  powers  of 
B  is  doomed  to  fail. 

This  failure  of  an  expansion  in  powers  of  B  was  demonstrated  by  an 
abortive  treatment.  The  eigenfunction  could  be  expanded  as  a  series 
of  Legendre  polynomials,  and  it  would  be  expected  that  the  coefficient  of 
Pj^(x)  would  be  the  largest.  The  fourth  derivative  of  a  Legendre 
polynomial  can  be  expressed  as  a  sum  of  Legendre  polynomials  of  orders 
up  to  k  -  4.  Hence,  equation  3.30a  can  be  written  as  a  set  of  equations 
for  the  coefficients  in  the  expansion,  and  the  boundary  conditions  yield 
two  further  relations  for  these  coefficients.  The  condition  that  this 
set  of  linear  homogeneous  equations  possess  a  solution  other  than  zero 
requires  that  the  determinant  of  the  set  vanish,  which  provides  an 
equation  for  the  eigenvalue  The  coefficients  of  the  expansion  then 

are  found  from  the  cofactors  of  the  determinant. 

This  is  a  reasonable  procedure,  but  it  is  found  that  the  terms  of 
the  determinant  increase  rapidly  away  from  the  diagonal.  Suppose  the 
determinant  is  truncated  at  some  value  N,  and  we  wish  to  find  the 
eigenvalues  to  first  order  in  B-  This  can  indeed  be  done,  with  a  result 
depending  on  N.  The  limit  as  N  approaches  infinity  can  be  taken,  and  the 
eigenvalue  does  approach  a  limit.  However,  if  we  attempt  to  calculate 
the  eigenfunctions,  also  to  first  order  in  B>  the  expansion  coefficients 
are  all  of  order  B  except  those  corresponding  to  N  -  1  and  N,  which  are 
of  order  unity.  Hence,  the  representation  clearly  has  no  well-defined 
convergence  process,  since  the  coefficients  may  be  made  to  change 
arbitrarily  by  changing  the  number  of  rows  kept  in  the  truncated 
determinant . 

We  have  developed  a  convergent  algorithm  for  computing  the 
eigenvalues  and  eigenfunctions  of  the  differential  equation  3,30a  with 
boundary  conditions  3.30b.  Before  embarking  on  this  analysis,  we  will 
prove  that  the  eigenvalues  of  3.30  are  all  real  and  positive.  This 
demonstrates  that  the  free  motion  is  stable,  that  is,  there  are  no 
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growing  oscillations.  The  proof  is  limited  to  small  displacements,  the 
range  of  validity  of  3.30,  but  that  is  the  only  range  of  interest. 

To  prove  the  eigenvalues  are  real,  suppose  w  is  a  complex 
eigenvalue,  with  eigenfunction  X.  Let  a  bar  above  a  symbol  denote  the 
complex  conjugate.  Then,  since  3.30  is  a  linear  system  with  real 
coefficients,  other  than  u  ,  we  see  that  w  is  also  an  eigenvalue,  and  X 
is  its  eigenfunction.  Multiply  the  differential  equation  for  X  by  X,  the 
differential  equation  for  X  by  X,  subtract,  and  integrate  over  x  from  -1 
to  1.  There  results 


-  ti>^)/dxXX  =  /dx 


■-jd  2.d  ^  pd^X 

X  — (1  -  X  )— X  -  P— 7- 
|dx  dx  ^^4 


dx 


3.32 


where  the  limits  are  implied.  Integrate  the  right  sides  by  parts.  For 
the  terms  independent  of  3,  a  single  integration  by  parts  yields  boundary 
terms  which  vanish  because  of  1  -  x  ,  and  integral  terms  which  cancel. 

For  the  P  terms,  two  integrations  yield  boundary  terms  which  vanish  by 
3.30b  and  its  conjugate,  and  also  integral  terms  which  cancel. 

Therefore,  the  right  side  of  3.32  vanishes.  The  integral  on  the  left  has 
a  positive  integrand  and  cannot  vanish,  so  the  only  possibility  is 
a  =  requiring  the  imaginary  part  to  vanish  and  u  be  real. 

To  prove  is  positive,  multiply  3.30a  by  Xj^  and  integrate  from  -1 
to  1.  Integrate  by  parts  as  above.  The  boundary  terms  will  vanish,  and 
there  will  result 


2  2  2 


/dx 


3.33 
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Since  all  integrals  are  positive,  is  positive.  The  same  procedure  can 
be  used  to  prove  the  solutions  are  orthogonal,  i.e.. 


/dxXj^X^  =  0  k  ^  I  . 


3.34 


We  shall  normalize  the  eigenfunction  so  that 


3.35 


These  results  are  all  derived  without  finding  the  explicit 
eigenvalues  and  eigenfunctions.  To  determine  the  solutions,  we  consider 
the  differential  equation  and  boimdary  conditions  with  the 
gravity-gradient  term  omitted.  Ve  call  this  the  bar  problem.  This  seems 
perverse,  since  we  are  dropping  a  term  with  unit  coefficient  and  keeping 
a  term  with  a  small  coefficient.  However,  this  simplified  problem 
contains  all  the  singular  properties  of  the  general  problem.  The 
eigenfunctions  of  the  bar  problem  form  a  basis  for  expansion  of  the 
complete  problem,  and,  whereas  we  may  expect  slow  convergence,  the 
solution  should  actually  converge,  unlike  the  previous  sltuation. 

The  bar  problem  has  been  considered  by  many  writers [13,  14,  and 
others] .  We  will  separate  the  eigenfunctions  into  the  two  classes  of 
even  and  odd  functions.  The  eigenvalues  of  the  bar  problem  will  be 
denoted  by  Vj^,  the  eigenfunctions  by  Yj^(x)  . 

The  bar  problem  is  defined  by  the  fourth-order  differential  equation 
and  boundary  conditions 


3.36a 
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2  3 

At  X  =  ±1  ^  =  0  .  3.36b 

dx"^ 


The  solution  of  3.36a  for  even  eigenfunctions  is 


Y,  =  B,  cosv,  X  +  C,  coshv,  x  . 
k  k  k  k  k 


To  satisfy  the  boundary  conditions,  we  require 


3.37 


3.38a 

3.38b 

These  equations  must  have  a  vanishing  determinant,  which  yields  the 
eigenvalue  equation: 


-Bj^cosVj^  +  Cj^coshVj^  =  0 


+Bj^sinVj^  +  Cj^sinhVj^  =  0 


tanv,  =  -tanhv,  . 
k  k 


3.39 


Tor  at  all  large,  the  right  side  tends  rapidly  to  -1,  so  Vj^  will  tend 
to  (k  -  l/4)ir.  We  have  the  explicit  sequence 


Vj/v  =  .7528094 

3.40a 

v^/if  =  1.7500053 

3.40b 

=  k  -  1/4 

to  8  decimals 

k  >  2  . 

3.40c 
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The  eigenfunctions  take  the  form 


=  Nj^[coshv^cosv^x  +  cosVj^coshVj^x  ] 


3.41 


where  denotes  a  normalizing  factor.  The  bar  eigenvalues  will  also  be 
normalized  to  unit  mean  square: 


/dxY,^  =  1  . 

k 


3.42 


The  integration  is  straightforward,  and  yields 


Nj^  “  C/2coshVj^|  cosVj^j  ) 


3.43 


For  the  odd  eigenvalues,  we  obtain  the  equation  analogous  to  3.39, 


tanVj^  =  tanhVj^ 


3.44 


with  the  solutions 


Vj^/TT  =  1.2498763  3.45a 

V2/'IT  =  2.2499998  3.45b 


k  +  1/4 


to  9  decimals 


k  >  2  . 


3.45c 
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The  normalized  odd  eigenfunctions  are 


\  =  N^lsinhvj^sinuj^x  + 


s  inhVj^x] 


3.46a 


=  l/(/2sinhVj^|  sinv^l  )  . 


3.46b 


let  us  expand  the  solutions  of  the  full  problem  3.30  in  terras  of  the 
solutions  of  the  bar  problem,  namely. 


oo 

Xj^(x)  =  . 


3.47 


The  boundary  conditions  are  automatically  satisfied.  The  differential 
equation  becomes 


Za,. 


3.48 


Multiply  by  Y  (x)  and  integrate.  The  left  side  reduces  to  the  term  in 

a  .  Define  a  matrix  element  coefficient  c.  by 
k,m 


9  ‘iY,  dY 

c  =  /dx(l  -  X  )- 


dx  dx 


3.49 
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Then  the  differential  equation  3.48  becomes  the  set  of  algebraic 
equations : 


+ 


2  2 
3“k 


Xc.  a,  =0  £,  =  1,. 

.  t.m  k.m 
m=l 


3.50 


where  the  prime  on  the  sum  means  the  term  m  =  £  is  to  be  omitted. 

The  evaluation  of  the  coefficient  c.  is  a  tedious  task  in 

jc  ^  in 

integration  of  elementary  functions.  Omitting  the  details,  we  have  for 
both  even  and  odd  functions 


c.  =  -32v  V  / (v  -  V  ) 
£,m  Jl  m  £  m 


I  ^ 


m 


3.51 


The  diagonal  elements  are  given  by 


2  2^  2 
=  3V^tan  V 


£ 


even 


3.52a 


'^£,£ 


2  2  ^2 


odd  . 


3.52b 


We  see  from  3.51  that  if  m  is  large  compared  to  £,  that  is,  far 

4 

above  the  diagonal,  the  coefficient  is  approximately  -32(£/m)  and  the 
terms  decrease  rapidly.  The  diagonal  elements,  which  are  the  largest, 
increase  as  £  ,  so  the  equations  may  be  normalized  by  dividing  by  the 
diagonal  elements  and  the  resulting  determinants  are  all  convergent. 

With  the  procedures  proven  legitimate,  we  did  not  attempt  to  treat 
the  eigenvalues  and  eigenfunctions  by  solving  the  determinants 
analytically,  but  resorted  to  a  computer  program  available  for  such 
purposes  at  Rand.  The  basic  program,  EIGRS,  was  prepared  and  is 
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supported  by  IMSL,  Inc.,  Houston,  Texas  (International  Mathematical 
Statistical  Libraries).  The  program  to  calculate  the  coefficients  ^ 
and  use  EIGRS  was  written  by  M.  D.  Lakatos. 

The  fundamental  modes  were  calculated  for  the  three  values  3=0, 

3  X  10~^,  lO'^.  The  infinite  set  of  equations  3,50  were  truncated  by 
multiples  of  5  up  to  40  equations,  and  the  eigenvalues  and  expansion 

coefficients  a^^  ^  were  determined.  The  first  20  even  and  odd  eigenvalues 
and  eigenfunctions  were  found  accurately  to  four  significant  figures  (no 
changes  of  10  ^  in  the  increase  from  35  to  40  equations),  and  the  lower 
ones  were  much  more  accurate  than  that.  The  first  four  odd  and  even 
eigenvalues  for  and  the  20th,  are  listed  in  Table  3  for  the  three 
values  of  3. 

We  may  immediately  draw  some  important  conclusions  from  Table  3. 

First,  in  spite  of  the  small  values  of  3,  even  the  fourth  eigenvalues 
experience  appreciable  changes  (5  percent) ,  and  the  20th  eigenvalues 
experience  40  percent  increases.  Second,  and  more  important,  the  lowest 
eigenvalue,  corresponding  to  the  first  even  mode,  is  very  nearly  an  integer, 
that  is,  the  lowest  vibrational  mode  is  at  nearly  three  times  the  orbital 
frequency.  If  the  perturbations  contain  terms  at  thrice  orbital,  and  it 


Table  3 


CHARACTERISTIC  FREQUENCIES 


3 

-5 

k 

0 

3  X  10 

10 

IE 

3.0000 

3.00032 

3.00107 

10 

4.2426 

4.2454 

4.2506 

2E 

5.4772 

5.4875 

5.5079 

20 

6.7083 

6.7380 

6.7916 

3E 

7.9376 

8.0057 

8.1157 

30 

9.1661 

9.2980 

9.4882 

4E 

10.3950 

10.6191 

10.9105 

40 

11.6248 

11.9709 

12.3808 

20E 

54.88 

64. 18 

76.82 

200 

56.37 

66. 16 

79.57 
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will  be  demonstrated  later  that  thermal  bending  effects  do  involve  such 
frequencies,  then  strong  resonant  responses  can  be  expected. 

A  computer  program,  also  written  by  M.  D.  Lakatos,  has  been  used  to 
calculate  the  eigenfunctions  from  the  expansion  coefficients  ^ 

and  the  bar  eigenfunctions  The  first  and  second  even 

eigenfunctions  are  plotted  in  Fig.  5,  the  fifth  even  eigenfunction  in 
Fig.  6.  Only  the  portion  for  positive  x  is  shown.  All  curves  are  for 
g  =  Also  graphed  on  Figs.  5  and  6  are  the  suitably  normalized 

Legendre  polynomials  P^,  and  In  Fig.  5,  the  departure  of  the 

true  eigenfunction  from  the  Legendre  polynomial  is  only  discernible  at 
the  extreme  right  end  of  the  curve,  that  is,  near  the  end  of  the  rod. 

The  Legendre  polynomials  have  strong  curvature  near  the  ends,  whereas  the 
true  eigenfunctions  are  straight  to  fourth-order  terras,  so  the  closeness 


Fig.  5  —  First  two  even  fundamental  modes 
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X 


j  Fig.  6  —  Fifth  even  fundamental  mode 


of  the  approximation  is  pleasing.  The  fifth  modes  show  signifioant 
departures  throughout,  reaching  a  difference  of  1.4  units  at  the  end.  As 
is  plausible,  more  of  the  energy  density  per  unit  length  of  the  modes  is 
toward  the  center  of  the  rod  for  the  true  eigenfunctions  than  for  the 
gr'avity“gradient  modes,  since  the  elasticity  suppresses  bending  at  the 
ends.  The  Nth  eigenfunction  has  N  zeros,  just  like  the  Legendre 
polynomial . 

With  the  eigenvalues  and  eigenfunctions  constructed,  we  are  now 
ready  to  determine  the  effects  of  initial  conditions  on  vibratory  motion. 
The  initial  values  of  the  coefficients  are  given  by 


A^(0)  =  fdxXj^Cxjw^Cx,  0) 


3.53a 
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A^(0)  = 

If  the  rod  is  started  from  rest,  A^CO)  =0.  We  will  return  to  this  point 
when  we  discuss  oiicrometeoroid  impacts,  which  impart  an  initial  velocity 
to  a  short  section  of  rod,  but  limit  ourselves  to  A^(0)  =  0  at  present, 
and  look  for  suitable  initial  conditions. 

Previous  investigations [2,  3,  11]  have  assumed  the  rod  is  distorted 
into  a  parabolic  shape.  However,  this  does  not  meet  the  boundary 
conditions,  and  the  ends  of  the  rod  would  jump  discontinuous ly  if  the 
constraints  establishing  the  initial  distortion  were  released.  We  choose 
as  our  first  initial  condition  the  simplest  even  polynomial  which  meets 
the  boundary  conditions,  has  vanishing  average  value,  and  is  normalized 
to  unit  deflection  between  center  and  ends.  This  proves  to  be  the  sixth- 
degree  expression: 


w,  =  (-29  +  105x 

O 


35x^  +  7x^)/77 


3.54 


2 

This  expression  and  the  comparison  parabola  (3x  -  l)/3  are  plotted  in 

Fig.  7.  The  two  curves  are  very  close,  the  maximum  difference  being  .043 
at  the  center  and  ends  (x  =  0  and  1),  and  .036  near  x  =  .65.  Hence, 
expression  3.54  is  a  proper  representation  for  a  near-parabolic 
displacement. 

Insert  3.54  into  3.53a,  and  replace  by  the  expansion  3.47. 

The  differential  equation  3.36a  may  be  used  to  replace  by 


and  the  resulting  expression  integrated  by  parts  seven  times.  The 
boundary  conditions  make  the  integrated  terms  vanish  for  the  first  six 
integrations.  Since  d  w^/dx  is  constant,  the  last  integration  evaluates 
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Fig.  7  —  Smooth  initial  conditions 


the  integral  as 


/dxWg(x)Y^(x) 


1.440 


llv. 


8  dx 


x=l 


3.55 


Evaluating  the  derivative  from  3.41  and  using  the  closeness  of  to 
(Jl  -  l/4)Tr  yields  for  3.55: 


1,440/2  (-1)^"^ 
11  7 


/dxw^(x)Y^(x)  = 


3.56 
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The  first  three  coefficients  have  the  numerical  values  .43955,  -.00122, 
.00005.  Thus,  almost  all  the  contribution  to  will  come  from  the  term 
in  the  expansion.  From  the  coefficients  as  determined  by  the 
computer  program,  we  have  ^  ~  .9966,  ^2  1  ~  ^3  1”*  The 

response  to  the  initial  displacement  3,54  is  therefore  given  to  better 
than  1  percent  accuracy  by 


W2(x,0)  =  .4381Xj^(x)cosw^9  -  . 0364X2 (xlcosw^B  - 


3.57 


The  motion  will  be  a  near-parabolic  shape  at  frequency  near  thrice 
orbital,  with  an  8  percent  admixture  of  double-noded  shape  (see  Fig.  5), 
at  frequency  near  5.5  times  orbital.  The  ’’throw  of  Xj^(x),  that  is,  the 
difference  between  maximum  and  minimum,  is  2.341,  so  the  throw  of  the 
first  term  is  1.0256.  We  observe  that  the  response  to  a  near  parabolic 
disturbance  is  primarily  an  oscillation  at  the  Initial  amplitude,  thrice 
orbital  frequency,  and  still  near -parabolic  shape. 

To  obtain  a  response  involving  many  modes,  we  use  a  "broken- line" 
initial  displacement.  This  is  equivalent  to  a  plucked  string.  The 
simplest  form  puts  the  break  at  the  center,  exciting  only  even  modes,  so 
we  choose 

i 

Wq(x)  =  .5  -  |x|  3--5i 


which  has  the  required  properties  (even,  break  at  center,  unit  amplitude, 
satisfies  the  boundary  conditions).  The  integral  involving  Wq  and  Y^(x) 
may  be  evaluated  easily,  yielding 


/dxwQ(x)Yj,  (x) 


2_  1 
^2  /2|cosVj^j 


COSVj^ 

coshv^ 


3.59 
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The  factor  multiplying  2/v^  is  1.1329  for  fc  —  1,  .9942  for  U  —  2,  and 
essentially  unity  thereafter.  We  see  that  the  terms  in  the  expansion 
with  respect  to  bar  eigenfunctions  of  A  (0) ,  the  initial  value  of  the 

K  2 

coefficient  of  the  true  eigenfunction,  only  fall  as  2/Vjj .  To  obtain  1 
percent  accuracy,  nine  terms  have  to  be  kept.  Thus,  the  coefficient 
Aj^(O)  is  approximately  given  by 


A,(0) 


2.2658 

2  ^k,l 

V, 


1.9884 


9 

+  2Za 


k,fc 


3.60 


With  the  coefficients  from  the  computer  program,  these  expressions, 
which  we  call  excitation  coefficients,  have  been  worked  out  using  the 
author's  HP-34C,  and  are  listed  in  Table  4. 

Again,  the  first  term  is  dominant,  but  the  higher  frequency  terms 
are  now  much  more  important  than  they  were  for  the  near-parabolic 
excitation  3.54.  The  shape  of  the  rod  at  the  initial  time  and  after  1 
and  10  orbits  has  been  calculated,  and  is  plotted  in  Fig.  8.  Because  the 
most  important  term,  that  at  nearly  thrice  orbital  frequency,  is  back  at 


Table  4 

EXCITATION  COEFFICIENTS  FOR  PLUCKED  ROD 


k 


1 

2 

3 

4 

5 

6 

7 

8 
9 


A^(0) 


.3954 

-.0886 

.0400 

-.0220 

.0147 

-.0096 

.0057 

-.0036 

.0031 
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Fig.  8  —  Response  to  initial  displacement  in-plane  i 


Its  initial  value  after  an  integral  number  of  orbits,  the  distorted  shape 
is  quite  similar  to  the  original  broken  line,  but  the  changes  are 
apparent.  Thus,  the  vibration  continues  at  the  original  amplitude,  but 
the  shape  will  now  be  a  complicated  function  of  time. 

The  analysis  has  all  been  for  in-plane  motion,  but  the  same 
techniques  adapt  directly  to  out-of-plane  motion,  with  equation  3.26 
providing  the  frequency  relation  and  everything  else  being  identical. 

The  investigation  shows  that  initial  tilting  or  distortion  will  yield 
libration  or  vibration  with  amplitude  substantially  equal  to  the  initial 
amplitude.  The  verticallty  and  straightness  requirements  must  therefore 
be  satisfied  by  the  initial  deployment  of  the  satellite.  The  verticallty 
should  not  be  difficult,  but  an  rms  lateral  displacement  of  .5  cm  over  a 
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length  of  1,500  m  promises  to  be  a  very  difficult  deployment  task,  since 
the  satellite  must  be  installed  in  a  dispenser  of  radius  limited  to  4  m 
by  the  payload  bay  dimensions  of  the  space  shuttle,  and  then  unrolled 
from  that  dispenser,  without  acquiring  any  permanent  bends  which  would 
lead  to  such  distortions. 

The  libration  and  vibration  would  decay  with  time  if  there  were  any 
damping.  It  has  been  shown  by  Burke [11]  that  the  "sticky  hinges"  of  the 
Stanford  design  will  be  locked  by  the  gravity-gradient  tension  and  will 
not  provide  the  required  damping.  Internal  friction  in  steel  is  also 
insufficient.  (The  damping  of  metal  strings  in  musical  instruments  is 
primarily  due  to  interaction  with  the  air.)  Viscous  damper  booms  like 
those  used  on  the  Radio  Astronomy  Explorer [5,  6]  cannot  be  attached  to 
PACSAT.  If  damping  is  to  be  provided,  some  other  procedure  must  be 
found. 

This  completes  our  discussion  of  the  basic  problem.  We  next 
consider  the  effects  of  various  perturbations,  beginning  with  orbital 
ellipticity. 
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IV.  ELLIPTICITY  EFFECTS 


There  have  been  previous  studies[7,  8]  of  the  effect  of  ellipticity 
on  the  libration  of  a  satellite.  We  shall  employ  a  different  technique, 
which  yields  the  same  driven  solution  and  stability  analysis,  and  is 
valid  for  both  in-plane  and  out-of-plane  motions.  We  shall  also 
investigate  the  effect  on  vibration,  which  will  prove  to  consist  of  small 
variations  of  frequencies  and  amplitudes. 

First,  let  us  consider  the  center-of-mass  motion.  This  will  be  an 
ellipse  in  the  equatorial  plane  with  the  earth  at  one  focus.  The  semi- 
major  axis  of  the  ellipse  is  denoted  by  a,  the  eccentricity  by  e,  and  the 
angular  coordinate  by  0,  measured  from  perigee.  The  eccentric  anomaly 
(see  any  book  on  celestial  mechanics)  is  denoted  by  E,  and  the  mean 
orbital  rate  by  n.  Then  we  have  the  equations: 


R  =  Ri_  =  R(i  COS0  +  i  sin9) 

R  X  y 

4 .  la 

R  =  a(l  -  e^)/(l  +  ecos0) 

4.1b 

COS0  =  (cosE  -  e)/(l  -  ecosE) 

4.1c 

nt  =  E  -  esinE 

4. Id 

2  3 

n  a  =  p  . 

4.1e 

Equation  4.1a  is  the  proper  generalization  of  3.1a;  4.1b  is  the 
equation  of  an  ellipse  of  the  indicated  dimensions.  The  angular 
coordinate  and  eccentric  anomaly  are  related  through  4.1c.  Kepler's 
equation  4. Id  expresses  the  time  in  terms  of  the  eccentric  anomaly,  and 
can  be  inverted  to  express  the  anomaly,  and  hence  the  angular  coordinate 
0,  in  terms  of  the  time.  Finally,  4.1e  is  the  generalization  of  3.5  to 
elliptic  motion. 
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The  libration  and  vibration  will  again  be  resolved  in  cylindrical 
coordinates.  While  equations  3.6a  and  3.6b,  which  give  the  first  time 
derivatives  of  the  unit  vectors  ij^  and  ig,  remain  valid,  equations  3.6c 
and  3.6d  must  be  generalized  to  include  the  variable  angular  velocity, 
hence, 


•2^ 

=  -e^i.. 


-->■ 


+  0i 


3  i 


3 

e  ■ 


The  libration  equation  2.33  becomes 


4.2a 

4.2b 


u 


1 


u 


2 


u 


3 


4,3a 


4.3b 


4 . 3c 


The  angles  ^  (in-plane  motion)  and  X  (out-of-plane  motion)  are 
introduced  by  3.4,  and  the  same  manipulation  used  to  separate  the 
equations.  The  angle  0  is  to  be  used  as  the  Independent  variable,  but 
the  transformation  is  somewhat  more  complicated.  It  is  found  that 


_M _  ^  1 

_3*2  1  +  ecos0 

R  U 


4.4a 


4.4b 
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9  _ 2esln9 

■2  1  +  ecos9 


4.4c 


The  equations  analogous  to  3.10  are 


1  ,2*2 

gPL  0 


..2 


cos^x(l  +  <P')^ 


1  +  ecos6 


2  2 

•(3cos  xcos  4/  -  1) 


4.5a 


if;" 


2esin6 


1  +  ecosO 


-if;’  + 


1  +  ecosO 


sini/icosif; 


2es  In  9 
1  +  ecos0 


+  2tanxx’  (1  +  il^') 


4.5b 


X  - 


2esln6 
1  +  ecDSO^ 


(1  +  ^.)2  +  3cos  If; 

VJ-  T  V  ^  1  +  ecosO 


sinxcosx  =  0  .  4.5c 


In  4.5a  0  is  regarded  as  variable.  The  most  important  difference 
between  4.5  and  3.10  is  the  inhomogeneous  term  on  the  right  of  4.5b. 
Although  equations  3.10  have  the  possible  solution  ^  =  x  =  0,  a 
permanently  vertical  rod,  this  is  not  possible  for  4.5.  This  driven 
solution  is  physically  very  reasonable.  The  velocity  vector  of  the 
center  of  mass  is  not  perpendicular  to  the  vertical,  so  an  in-plane 
attitude  oscillation  can  be  expected. 

The  intention  is  to  launch  PACSAT  into  a  circular  orbit,  so  it  can 
be  expected  that  e  will  be  small.  The  previously  cited  Radio  Astronomy 
Explorer  satellite  (1968-55A),  whose  orbit  and  physical  character  is  most 
similar  to  PACSAT,  had  an  orbital  eccentricity  of  .0004.  An  eccentricity 
of  .01  at  an  altitude  of  one  earth  radius  corresponds  to  a  velocity 
injection  error  of  183  fps,  far  above  any  reasonable  value.  We  will 
calculate  the  amplitude  and  frequency  of  the  motion  to  second  order  in  e, 
which  should  be  more  accurate  than  necessary.  The  amplitude  can  be  found 
by  expanding  4.5b  and  4.5c  to  second  order,  but  the  frequency  requires 
some  third-order  terms.  The  initial  conditions  will  be  left  unspecified, 
since  PACSAT  may  be  injected  at  an  arbitrary  point  along  the  orbit.  To 
second  order,  equations  4.5b  and  4,5c  become 
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ti 


2esine  ,  3  _  2esin9  , 

1  +  ecDso'^  1  +  ecos9*^  1  +  ecos9 


4 . 6a 


X 


fr 


2esln9  , 
1  +  ecosO^ 


1  + 


3 

1  +  ecos9 


4.6b 


The  first -derivative  terms  may  be  eliminated  by  the  transformation 


<l>  =  f  /  (1  +  ecosG) 


4.7a 


X  =  g/(l  +  ecosG) 


4.7b 


which  yields  the  simpler  equations,  correct  to  second  order. 


£"  +  (3  -  2ecos0)f  =  2esin0  +  2gg'  4.8a 

g"  +  (4  -  2ecos8)g  =  -  2gf'  .  4.8b 


Separating  the  first-order  and  second-order  parts  yields 


f  =  + 

^2 

4.9a 

g  =  gj  + 

4.9b 

q  +  3f^ 

=  2esin6 

4 . 9c 

si  +  4s^ 

^  0 

4.9d 

±2  +  3f2  =  +2ecos0f^  +  2g^g| 


4.9e 
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“  +2ecoseg^  -  ' 


4.9f 


The  solutions  of  4.9c  and  4.9d,  subject  to  the  initial  conditions,  are 


=  esin0  +  A^cos/36  +  Bj^sinv^B 


4.10a 


®1 


C^cos28  +  Dj^sin20  . 


4.10b 


Thus,  the  driven  part  of  the  in-plane  motion  has  amplitude  e  and 

oscillates  at  the  orbital  frequency,  whereas  the  initial  conditions 

provide  an  oscillation  at  the  natural  frequency  'fb  times  orbital  with  the 
2  2  1^2 

amplitude  (A^  +  B^)  '  .  There  is  no  driven  out-of-plane  motion,  only  the 
natural  double-frequency  oscillation. 

The  right  side  of  4.9e  Involves  terms  in  (*^  -  1)0,  +  1)0,  20, 

and  40.  The  solution  is 


^2  "  iieAj^[(2/3  +  l)cos/3  -  1)6  -  (2/3  -  l)cos(/3  +  1)6] 
+  YxeB^{{2/3  +  l)sin(/3  -  1)9  -  (2/3  -  l)sin(/3  +  1)9] 

+  if  H  -  ifsV”*® 

2 

-  e  sin26  +  A2COS/J6  +  B2sin/36  . 


4.11 
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The  right  side  of  4.9f  involves  terms  in  (2  +  /3)5  and  (2  -  v^)0,  with 
solution: 


®2  13 


(4  -  /3)(A^D^  +  B^C^)cos(2  +  /3) 


+  (4  +  /3)(A^D^  -  B^C^)cos(2  -  i^)e 
-  (4  -  ’/3)CAj^Cj^  -  Bj^D^)sin(2  +  /3)0 


-  (4  +  /3) (A^C^  +  B^D^)sln(2  -  /3)e  +  C2Cos2e  +  D2Sin2e 


.  4 


To  second  order,  the  out-of-plane  motion  is  independent  of  the 
eccentricity  except  as  e  appears  in  the  constants  A^,  which 

involve  the  initial  conditions. 

Thus,  the  amplitudes  of  the  in-plane  and  out-of-plane  motions  have 
been  obtained  to  second  order.  To  determine  the  frequencies,  it  is 
necessary  to  expand  4.5  to  third  order.  The  expressions  for  f2>  gj^> 

and  are  substituted,  and  the  terms  of  apparent  frequency  /3  and  2  are 
collected.  The  combined  amplitude  of  these  terms  provide  the  correction 
to  the  frequency.  Thus,  indicating  only  the  third-order  parts  on  the 
right,  we  have 


f"  +  (3  +  e^)f  =  2ecos6f2  "  e^cos20fj^ 


+  2[gj^g^  +  g^g^  -  ecosBg^gJ  +  esinOg^  +  +  fj]  4.13a 


.12 
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g"  +  (4  +  e‘^)g  =  2ecos9g  -  e  cos20g 


-  2[g^f2  +  S2^i  ~  ®cos0g^f|  +  eslnef^g^ 
^  lx. 2  3^2  4  3, 

+  gj  -  • 


4.13b 


Terms 


in  Aj^cos  /30 


collect  on  the  right  side  of  4.13a  as 


35  2 

U® 


3  2 

+  |(A^ 


+  B^) 


13''^!  ^ 


4.14 


from  which  we  deduce  the  modified  in-plane  libration  frequency 


2  _  ,  24  2  _  3,, 2  2  2  2. 

“iP  ^  "  TT^  "  13^^1  ^ 


4.15a 


IP 


=  /3 


,  4  2  1,.2  ,  „2. 

^  11®  4^^1  ®1^ 


12  2 
I3<'^1  ^  °1^ 


4.15b 


Similarly,  collecting  terms  in  cos20  on  the  right  side  of  4.13b  yields 


+  j~(A^  +  B^)  -f-  +  D^) 


4.16 


from  which  we  deduce  the  modified  out-of -plane  libration  frequency 


2  /  6,. 2  ^  2.  18  2  2. 
=  -  13^^  +  ®1^  "  13^S  + 


“op  =  2 


1  -  ^(a2  b2)  -  ^icl  4-  d2) 


4.17a 


4.17b 
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The  expressions  4.15b  and  4.17b  are  generalizations  of  the  results 
of  Refs.  7  and  8,  which  treat  only  the  in-plane  motion.  We  see  that  the 
frequency  of  the  out -of -plane  motion  is  independent  of  e  to  second  order, 
except  through  the  Implicit  dependence  via  the  initial  conditions. 

The  simplest  case  which  retains  full  physical  meaning  is  for  the  rod 
to  start  at  rest  at  perigee  with  an  Infinitesimal  displacement  in  an 
arbitrary  direction.  Then  the  constants  become 


The  frequencies  become 
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X  =  XoCosUgpO  + 


^o|” 


—cos  (2  +  »^)e 

/3 


+  (4  +  /3) 


-^cos(2  -  /3)6  -  (J^«siti(2  -  /3)6 
/3  ° 


+  5iJjQsinu)|^pB  >  -  -eXpCcosO  +  cos  39)  . 


+  iJ;Qsin(2  +  /3)6 
+  llecosw^pB 

4.20b 


For  most  reasonable  values  of  eccentricity  and  initial  displacement, 
the  second-order  terms  are  very  small.  The  tension  is  given  to  first 
order  by 


Tg  =  |pL^0^[l  +  2/3(-ApSin/3e  +  BpCOs/To)] 


4.21 


where  the  e  dependence  is  in  0  . 

There  is  a  question  concerning  the  stability  of  the  out-of-plane 

3 

libration.  Let  equation  4.5b  be  linearized  by  dropping  the  X  terms, 
initial  deflection  of  1  meter  at  the  tip  of  the  rod  corresponds  to 
X  =  .00133,  so  expressions  of  the  order  of  10  ^  are  being  neglected. 
Further,  to  simplify  the  analysis  while  retaining  the  significant 
effects,  let  the  initial  conditions  be  chosen  so  =  0  (start  the 

in-plane  motion  from  rest,  90  deg  from  perigee,  with  initial  deflection 
e)  .  Then  we  have  the  simplified  form  for  jfi: 


An 


3  2 

esin9  -  sin29 


4.22 


and  the  differential  equation  for  g  becomes,  to  second  order  in  e. 


g"  +  [4  -  3e  cos20]g  =  0 


4.23 
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Equation  4.23  is  in  the  standard  form  of  the  Mathieu  equation[ 17 ] : 


y"  +  (a  -  2qcos20)y  =  0  , 


4.24 


The  stability  theory  of  such  equations  is  developed  in  Ref.  17,  Chapter 

4.  It  raay  be  summarized  as  follows.  When  q  =  0,  equation  4.24  will  have 

periodic  solutions,  with  period  Ti,  if  a  is  the  square  of  an  even  integer. 

If  q  is  not  zero,  there  will  be  two  specific  values  of  a,  regarded  as  a 

function  of  q,  such  that  the  equation  has  a  solution  of  period  ir.  The 

two  values  correspond  to  even  and  odd  solutions.  These  values  of  a,  call 

them  a  and  a„,  then  describe  two  curves  in  the  a,q  plane,  which  issue 
E  O'  2 

from  the  a  axis  at  a  =  4k  ,  k  integral.  If  the  value  of  a  which  actually 
appears  in  the  equation  4.24  lies  between  the  pair  of  curves,  the 
solution  of  4.24  is  unstable;  if  a  lies  outside  the  region  between  the 
curves  the  solution  is  stable.  According  to  Ref.  17,  p.  17,  the  curves 
Issuing  from  a  =  4  are 


Ug  =  4  +  +  O(q^)  4.25a 

Uq  =  4  -  y|q^  +  OCq"^)  .  4.25b 


We  see  that  for  4.23  the  value  of  a,  4,  is  at  the  junction  of  the  curves, 
and  the  problem  of  stability  is  not  resolved  to  the  second  order  and 
should  be  carried  to  the  fourth  order. 

If  we  assume  X  is  infinitesimal,  the  solution  of  4.5a  is  found  to 


fourth  order  as 


4.26 


3  2  3/13  3 

il;  =  esinS  -  -^e  sln29  +  e  '^sinG 

+  e^|-^sin29  - 


and  when  this  is  substituted  into  4.5b  and  the  lengthy  manipulations  are 
carried  out,  the  fourth-order  generalization  of  4.23  is 


4  -  -  3e^cos26  +  ye^(cosO  —  cos30) 

o  4 


.  4/33  .  1,143 

+  e  l-^cos20  +  ■  ■^q^cos401 


g  =  0  . 


4.27 


This  equation  is  in  the  form  of  a  generalization  of  Mathieu's 
equation  known  as  Hill's  equation  (Ref.  17,  Chapter  6).  The  stability 
theory  proceeds  in  the  same  manner.  The  presence  of  the  odd  harmonics  of 
9  requires  that  the  periodic  solutions  have  period  2if.  After  extensive 
calculations ,  the  curves  a^,  and  are  found  to  be 


a„  = 


,  237  4 


^  13® 


4.28a 


4.28b 


and  the  solution  of  4.27  is  unstable.  However,  if  e  —  .01,  it  will  take 
about  10^  orbits  for  the  instability  to  produce  observable  effects,  and 
this  is  over  4,500  years.  Hence,  although  the  out-of-plane  motion  is 
nominally  unstable,  there  will  be  no  effects  during  the  operating 
lifetime  of  PACSAT. 
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We  next  consider  the  effect  of  ellipticity  on  vibration.  The 
equations  3.20  must  be  generalized.  We  will  keep  only  first-order 
changes  produced  by  the  ellipticity.  The  differential  equation  2.34, 
expressed  in  cylindrical  coordinates,  becomes 


4.29a 


4.29b 


4.29c 


If  we  multiply  the  first  equation  by  u^,  the  second  by  the  third 

by  u^,  and  add,  we  obtain  an  expression  for  the  perturbation  in  the 
tension.  After  extensive  simplification,  it  becomes  to  second  order 


■  ■  ■  2 

1 - 1  =  _2(6  +  +  xw,)  -  2xw  +  29  . 

9s  p  ^  ^  ^ 


4.30 
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The  important  difference  between  4.3  and  4.29  is  that  in  4.3b,  the  term 
’Bu^  is  first  order  and  contributes  the  driving  term,  but  in  4.29b  the 
term  'Bw^  is  third  order,  since  is  second  order.  Hence,  there  is  no 
driving  term  in  4.29.  In  4.29b  and  4.29c,  the  expression  4.30  can  be 
substituted  in  the  first  term  on  the  right,  and  we  see  that  it 
contributes  only  second-order  quantities,  (In  all  this  discussion,  e,  </i, 
X,  w^j  Wg  are  regarded  as  first  order.)  We  define  the  convenient 
quantities : 


F  =  w^/Cl  +  ecosB) 


4.31a 


G  =  Wg/ (1  +  ecosB)  . 


4.31b 


We  are  not  interested  in  the  interrelation  between  the  ellipticity 
and  the  initial  conditions  on  the  librations  as  they  affect  the 
vibration.  Hence,  we  will  choose  the  initial  conditions  so  that  ~ 

—  Q  —  =  0  (see  equations  4.10)  .  Then  the  equations  for  the 

vibration,  complete  to  first-order  corrections  due  to  ellipticity,  become 
analogous  to  3.27: 


|(G''  +  G)  -  |j(l  -  -  6(1  -  4«058)^ 

with  the  same  boundary  conditions  as  in  3.27.  If  the  first-order  and 
second-order  parts  of  4.32a  are  separated  (analysis  of  4.32b  is 
superfluous  for  the  same  reason  as  in  3.20),  the  equations  for  the  first 
order  part  reduce  to  3.27.  The  second-order  part  is  the  forced 
oscillation  governed  by 


4.32a 

4.32b 
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4 

8  F. 


2^.,  8  ,  o  2  ,  a  Q  1 

— F"  -  —(1  -  X  )- - V  f5 - 7  =  4e3cose 

3  2  3x'  ^3x  ^4 


3x 


3x 


4  • 


4.33 


Substituting  for  the  expansion  3.28  and  for  a  corresponding 
expansion  with  replacing  Aj^,  there  results: 


CO  ^  3^ 

z|(b;^  +  =  4e6cose2A^(e)— ^ 

B"  +  w.^B,  =  6e3cos0ZA.  (0)cosci3  el  dxX.  (x)— r  . 

k  k  k  ^  fc  ^  J-1  3x^ 


4.34 


4.35 


The  integral  on  the  right  side  of  4.35  may  be  expressed  in  terms  of  the 
expansion  coefficients  of  the  true  eigenfunctions  with  respect  to  the  bar 
eigenfunctions  (see  equation  3.47)  as 


Za,  a„  V 
^  k,m  t  ,m  la 


4.36 


from  which  we  obtain  the  solution  for 


B 


k 


{» 


cosCm^ 


cos +  1)6 

-  (O,  +  1)^ 


4.37 


There  are  no  possible  vanishing  denominators  in  4.37,  so  there  are 
no  resonances,  and  the  perturbations  remain  small.  We  see  that  the 
effect  of  ellipticity  on  vibration  is  to  provide  only  small  corrections. 
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The  changes  in  the  vibration  frequencies  are  of  order  e  and  may  be 
neglected  completely. 

This  completes  the  analysis  of  ellipticity  effects.  We  next 
consider  effects  of  external  forces,  beginning  with  the  oblateness  of  the 
earth . 
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V.  EARTH  OBLATENESS  EFFECTS 


The  analysis  of  the  earth  oblateness  is  quite  complicated. 

Oblateness  affects  the  libration  and  vibration  equations  in  two  ways: 
the  direct  effect  caused  by  the  VF  term  in  equations  2.33  and  2.34,  and 
the  indirect  effect  caused  by  the  change  in  the  orbit  of  the  center  of 
mass.  Furthermore,  it  develops  that  there  is  a  resonance  in  the  out- 
of-plane  libration,  so  many  terms  must  be  kept  in  coefficients  to  make 
sure  that  all  significant  terms  are  included. 

The  investigation  begins  by  finding  the  modified  orbit  of  the  center 
of  mass.  As  is  well  known,  oblateness  causes  the  orbit  plane  to  process 
around  the  pole,  and  causes  the  perigee  to  process  aroixnd  the  orbit.  The 
equations  of  libration  and  vibration  must  be  expressed  in  terms  of  axes 
moving  with  the  precessing  plane  if  the  motion  is  to  be  understood.  They 
are  expressed  thusly,  and  the  resonance  treated  carefully.  The  vibration 
analysis  will  be  very  simple,  since  there  are  no  resonances  to  the  order 
treated . 

The  analysis  of  the  center-of -mass  motion  will  follow  Brouwer  and 
Clemence[ 18 ] .  It  would  have  been  much  easier  if  we  could  just  write  down 
the  results.  However,  although  the  precession  frequencies  are  matters  of 
general  knowledge,  we  have  not  found  a  reference  that  states  the 
amplitudes  of  the  oscillatory  terras  to  the  desired  accuracy.  We  have  had 
to  derive  them,  and  do  not  mean  to  imply  that  we  have  used  the  simplest 
procedure. 

Consider  a  satellite  moving  in  a  near-circular,  near-equatorial 
orbit  of  eccentricity  e  and  inclination  I.  Let  the  latitude  of  the 
satellite  be  X,  the  longitude  tfi.  We  express  the  coordinates  in  terms  of 
Delaunay  variables  (Ref.  18,  p.  541).  The  mean  anomaly  (nt  in  equation 
4. Id)  is  denoted  by  i,  the  argument  of  perigee  by  g,  and  the  longitude  of 
the  node  by  h.  Then  the  latitude  and  longitude  are  given  by 


sinX  =  sinlsin(0  +  g) 


5.1a 


-  68 


=  h  +  tan  [cosltaii(0  +  g)  ] 


5.1b 


The  radius  R  and  angular  coordinate  0  are  still  given  by  4.1b-e. 

The  coordinates  will  be  expanded  to  third-order  terms  in  the 
eccentricity  and  inclination.  The  result  is 


R  =  a 


1  + 


12  /  J  Jv 

2®  "  8®  ) 


3  3\ 


12  3  3 

cosil  -  -xe  cos2i  -  cos3Jl 
L  o 


5.2a 


sin£  =  I 


|l  -  e^  -  -^I^IsinC?.  4-  g)  +  e[sln(2A  +  g)  -  sing] 


+  |e^sin(3t  +  g)  +  -le^sinCA  -  g) 


5.2b 


=  £  +  g  +  h  +  |2e  -  4-  ■^e^sin2£.  +  ■— e 


5  2 


13  3 


sin3£ 


-  ^I^[sin(2£  2g)  +  2esln(3£  +  2g)  -  2esin(£  +  2g)]  .  5.2c 


If  there  were  no  perturbations,  i  would  vary  linearly  with  time,  and  all 
the  other  elements  would  be  constant.  The  disturbing  force  induces 
element  variations,  which  will  be  treated  as  first-order  quantities 
because  they  are  proportional  to  at  least  the  first  power  of  the  strength 
of  the  disturbing  force.  Denoting  the  variations  by  6a,  6e,  etc.  ,  we 
have  for  the  variations  of  the  coordinates ; 


69 


-6R  =  — —  +  e-(l-  ^e^\cosJl  -  ecos2Jl  -  ■^e'^COsSi!,  6e 

a  a  a  \  8  I  o 


9  2 


+ 


(1  -  1^^) 


+  esinZJ.  +  ■^e^sin3£  eS£ 


5.3a 


6(sin£)  =  -—sinj:  +  I 


-2esin(£  +  g)  +  sln(2£  +  g)  -  sing 


+  -^eslnCS?.  +  g)  +  -^esinC^  -  g)  6e 
+  l|cos(£  +  g)  +  e[cDs(2£  +  g)  -  co3g]|5(£  +  g) 

+  l|cos(2£  +  g)  +  cosg  +  |ecos(3£  +  g)  +  ■|ecos(£  -  g) 


ei5£ 


-  -l-I^sinCft  +  g)(SI 


5.3b 


6<}>  =  6(£.  +  g  +  h)  + 
1  2 


2  ^1  -  -le^lsinil  + -|esin2£  +  sin3£  6e 


+ 


5  13  2 

cos£  +  •::ecos2£  H - re  cos3£ 

2  4 


e6£  -  il6rsln(2£  +  g) 


-  yl^|[sin(3£  +  2g)  -  sin(£  +  2g)]6e  +  cos(2g  +  2g)6{£  +  g) 

^  I 


+  [cos(3£  +  2g)  +  cos(£  +  2g)le6£ 


5.3c 


The  terras  of  5.3  have  been  arranged  to  eliminate  apparent  singularities 
in  the  perturbations  of  the  elements  when  e  and  I  are  small. 

The  disturbing  force  is  the  gradient  of  a  potential: 

2 

^■^2^e  2 

- ^  (1  -  3sin  £) 

2R 


U  = 


5 . 4a 


5.4b 


f  =  VU 


2R^ 


2  ->■ 
[(1  -  3sln  £)ij^  +  2sini;cos£ijj)] 


-3 

The  oblateness  coefficient  has  the  numerical  value  1.0826  x  10  ,  and 

R  is  the  radius  of  the  earth.  The  theory  as  developed  in  Ref.  18  uses  a 
e 

Hamiltonian  formalism.  Coordinates  L,  G,  H,  canonically  conjugate  to  Z, 
g,  h,  are  introduced,  and  the  Hamiltonian  is  represented  in  terms  of 
them.  A  canonical  transformation  is  performed,  which  introduces  a  new 
set  of  elements  such  that  the  dependence  of  the  Hamiltonian  on  g,  h 
disappears  to  first  order  in  J^.  The  variations  of  the  elements  are 
given  in  a  closed  but  very  clumsy  form  on  p.  568  of  Ref.  18.  Expanding 
to  second  order  in  e  and  I  yields  the  variations  in  the  elements  as 
required  by  equation  5.3,  as  follows: 


=  |.j[4ecosil  +  6e^cos2i.  +  2I^cos(2A  +  2g)  } 


5.5a 


2  -  yl^jcosS.  +  3ecos2£  +  -^e^co 
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3. 

Se  =  p. 


'+  cos(£  +  2g)  +  ^cosC3£  +  2g) 


s3)l 


5.5b 


3  7 

61  =  ^jl  cos(2Z  +  2g)  -  ecos(£  +  2g)  4-  •^ecos(3&  +  2g) 


5.5c 


2|l  -  -  --I^jsinX  +  3esln2A  +  -j|e^sin3i!. 


=  -|a< 


-  V 
2 


sin(J!,  4-  2g)  -  -^sinOil  +  2g) 


5.5d 


1(6)1  +  6g)  =  Ijl 


13esin)l  -  sin{2J,  +  2g)  +  esln()l  +  2g) 
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6  (£  +  g  +  h)  =  -^jiyesinJl  +  6e^sln2£  +  r^s±n(2ii,  +  2g)]  5.5f 

j  =  J^Rg/a^  =2.7  X  10  ^  [PACSAT  at  1  earth  radius]  .  5.5g 


Where  these  are  substituted  into  5.3,  there  is  very  extensive 
simplification.  The  expressions  which  finally  result  for  the  coordinates 
themselves  are  given  to  second  order  by 


T  3.  ,  1  2/,  9.\  ,  9.^2  /t  3.\ 

1  +  2^  (1  -  2J)  +  -  e(l  -  ^.l) 


3.1 


cosl 


R  =  a7 


5 . 6a 


-  |•j)cos2)^  +  -131^008(211  +  2g) 
slni;  =  I  1^1  -  pjsinCJl  +  g)  +  e(l  +  ^jjsin(2Z  +  g)  -  e|l  +  Yjlsingj 


5.6b 


(|)  =  £.  +  g-t-h  +  2e|l+^j  Jslnt  +  +  ~2^s±rL2l 


-  +  j)I^sin(2il  +  2g)  . 


5.6c 


The  perturbed  frequencies  are  given  on  p.  572  of  Ref.  18  as 


3. 

/,  3  2 

3^2\ 

1  +  -^ 

(l  +-e 

2^  ) 

g  =  3nj (l  +  26^  - 


5.7b 
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h  = 


3  /  2 

-2^3(1  + 


5.7c 


2 

n  a  =  li  . 


5.7d 


The  frequencies  may  be  interpreted  as:  H  is  the  perigee  to  perigee 

•  •  *  •  •  . 

or  anomalistic  frequency,  )l  +  g  is  the  nodal  frequency,  and  £  +  g  +  h  is 

the  longitudinal  recurrence  frequency.  The  changes  indicated  by  5.7  show 
that  the  perigee  advances  and  the  node  regresses.  It  is  convenient  to 
redefine  the  orbital  elements  a,  e,  and  1,  which  are  arbitrary  constants, 
to  simplify  5.6.  We  shall  choose  a  new  set  of  parameters  such  that  the 
constant  term  and  the  term  in  cost  in  5.6a,  and  the  term  in  sin()l  +  g)  in 

5.6b,  have  the  same  coefficients  as  in  5.2a  and  5.2b.  With  the  new 

parameters,  we  obtain 


R  ==  ajl  +  -  ecosS.  -  ^1  -  -^jcos2J!,  +  •^I^cos(2Jl  +  2g) 


1.^2 


sin£  =  I 


sin(£  +  g)  +  e|l  +  |-j)sin(2£  +  g)  -  e  +  -|jjsing 


5.8a 


5.8b 


=  £  +  g  +  h  +  2e|l-*-'2j) 

-  +  |a  )sin(2£  +  g) 


sin£  +  e' 


!/5  ,  13. \ 

(4  "-“8^) 


sin2£. 


£  =  n 


5.8c 


5.8d 


i  +  g  =  nfl  +  13(3  +  5e^  - 


5.8e 
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S.  +  g  +  h  =  n 


1  + 


5.8f 


The  angular  rate  in  orbit  0  is  given  by 


0  =  n 


^  13(1  +  +  2e^l  +  -|•jjcosi  +  e^(-|  +  ^jjcos2£ 


3.^ 


-  |-JI^cos(2fc  +  2g)  . 


5.9 


With  the  results  5.8  and  5.9,  which  describe  the  motion  of  the 
center  of  mass  in  the  presence  of  oblateness,  we  can  proceed  to  determine 
the  libration  and  vibration  equations.  We  first  must  introduce  a 
coordinate  system  to  represent  the  in-plane  and  out-of-plane  motions. 

This  will  be  chosen  as  the  radial  unit  vector,  the  local  unit  horizontal 
velocity  vector,  and  the  cross  product  of  the  first  two.  Expressing 
these  coordinates  in  terms  of  the  latitude  and  longitude,  we  obtain 


+  cosJC^i  )/0  5.10b 

i„  =  (cosX^i.  -  Xi  )/6  .  5.10c 

J  <L  0 


The  displacements  u,  w  will  be  resolved  into  components  in  the 


directions  defined  by  5.10,  thus 
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u  = 

C11U. 

w  =  Wiii  +  W2i2  +  V3  • 

The  terms  Wj  are  the  radial  components,  u^,  are  the  in-plane 

components,  and  n^,  are  the  out-of-plane  components.  As  before,  is 
close  to  unity,  is  second  order,  and  the  others  are  first  order. 

The  time  derivatives  of  the  unit  vectors  are  needed  to  express  the 
acceleration  in  components.  We  introduce  a  quantity  P,  which  we  will 
call  the  twist,  by  the  relation 


P  = 


*2  6  1  ' 

cosJi!sin£(|)  -  —  /$cosJC  . 

6  J 


5.12 


The  time  derivatives  of  the  unit  vectors  are  given  by 


i 


1 


5.13a 


-0il  +  Pi3 


5.13b 


5.13c 


We  see  that  the  twist  P  measures  the  internal  rotation  of  the 
coordinate  system.  It  is  easily  shown  from  the  equations  of  motion  that 
P  vanishes  for  a  purely  elliptic  orbit.  For  motion  under  oblateness 
forces,  the  various  derivatives  in  5.12  may  be  reduced  using  5.8b  and 
5.8c,  and  the  result  to  second  order  is 


P  =  -3njl  |^sin(7.  +  g)  +  |esin(2i  +  g)  +  ^esingj  . 


5.14 
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The  acceleration  and  gravity-gradient  terras  combine  to  form  the  three 
expressions  as  left  sides  of  equations  analogous  to  4.3, 


5.15a 


5.15b 


5.15c 


and  similar  expressions  for  w.  Since  P  is  proportional  to  jl,  which  is 
-5  2 

on  the  scale  of  10  ,  the  terms  in  P  will  be  neglected.  The  term  PBu^^ 

in  5.15c  shows  that  the  twist  induces  a  driving  term  in  the  out-of-plane 
libration. 

We  next  have  to  transform  the  dyadic  VF  •  u  on  the  right  side  of 
2.33.  We  observe  from  5.4b  that  the  force  vector  F,  evaluated  at  an 
arbitrary  point,  r,  JJ,  0,  can  be  written  in  the  form: 


F(0  =  F^Cr.JOi^  +  . 


5.16 


The  displacement  from  the  center  of  mass  can  be  written  in  terms  of  the 
spherical  coordinates,  and  also  in  tha  form  5.11.  Using  V  as  the  generic 
term  for  the  displacement,  we  have 


V  =  V,  5,17a 

r  1 

V|,  =  (iv^  +  cosW^l/e  5.17b 

=  (cosj;0V„  -  iv,)/9  . 

9  z  J 


5.17c 
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The  changes  in  latitude  and  longitude  are  given  by 

5X  =  V  /R  5.18a 

6oi  =  V^/Rcos£  S.lSb 


and  the  variations  in  the  unit  vectors  i^. 


by 


^^r  =  (Vi; 

5.19a 

"  "(Vl 

5.19b 

=  '('"A  ■  • 

5.19c 

In  equations  5.17  -  5.19,  the  quantities  Z,  R,  and  the  time  derivatives 
are  evaluated  at  the  center  of  mass.  We  now  have  for  the  variation  of  F, 


5F  = 


3F  dF. 


1  +  F,6i  + 

r  1  r 


5.20 


After  extensive  reduction,  this  leads  to  the  form: 
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6F  = 


+ 


3F  /9F 

+  E\3x'  -  ''2/“''2  + 


\Q^  -  "2)“' 

S2  I 


.2  3Ft  ■  r.  (.12 
1  sinXcosXq) 


■^2^2 


/3F 

+  ^  +  tan£FjV3 


+ 


R0 

•  3F  *  •  /3F 

co5^ 2  — 2  tani!F  )V 

3R  1  „i2  \3£  tani^a^jv^ 


0 


R0 


ifr  -t-^ 
^Vi  e' 


•2  3F 


3£ 


■(ax  tani!F2)’' 

•2  \ 

)v  1 
;2  2/  3  3 


5.21 


We  will  keep  only  first-order  terras  in  e  and  I.  On  substituting  from 
5 . 4b  and  5.8,  we  have 


yj  8lsin(&  +  g)V3]i3^  -  [V2  +  2Icos(Jl  +  g)V3]i2 

af  =  1  !!■  2_e 

2  „5 


R 


[3V3  -  8rsin(Jl  +  g)V^  +  21cos{«.  +  g)V2]i3 


5.22 


If  we  keep  the  same  order  of  accuracy  in  the  expressions  5.15,  we  obtain 
the  equations  of  libration; 


u^  -  n^  jd  +  +  10e(l 

+  2en^(l  -  ■|:jjsln^.U2 
-  n^  -  ^^jecosJlu^ 

+  bnjlsinCi.  +  g)^^  + 


+  3j)cost 


-  2n(l  +  |j  +  ecos£ ju2 
8T, 


2  “0 
-  15n  jlsInCt  +  g)^^  =  ~ — 5.23a 

pL 

+  2n  +  Zecos^ju^  -  2en^|l  -  ^jjsin£Uj 

2 

6n  jIcos(2.  +  g)^3  =  — ^2 

pL 


5.23b 
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+  n 


1  +  4-  3  jecosA ju^  -  15n^jIsin(Jl  +  g)^2^ 


2  ^'^0 

-  6njlsin(i.  +  g)yi^  -  3n  jIcos(£  +  ^  ~i“3  * 


5.23c 


We  apply  the  same  transformations  to  these  equations  as  used  in 
3.10,  except  that  the  time  is  retained  as  the  independent  variable. 
Keeping  only  first-order  terms,  we  find: 


5.24a 


Tq  =  •^L^|n^|3  +  +  10e{l  +  3j)cos)i  +  2n 

-  2en^  -  ■|^jjsin^i/'  +  15n^jlsin(<^  +  g)x^ 

+  n  js  |l  +  -I-  (11  +  54j )ecos^^ I 'J' 

=  2en^  |l  -  •^■j|sinJl  -  6njIsin(Ji  +  g)x  -  6n^jIcos{^  +  g)x  5.24b 

X  +  n^[4{l  +  3j)  +  (15  +  69j)ecos£.]x  ■=  15n^jIsln(A  +  g) 


2  •  2 
+  15n  jelsin(2£  +  g)  +  6njlsin(£  +  g):|^  +  3n  jIcos()l  +  g)i*t  .  5. 


We  see  that  the  in-plane  and  out-of-plane  frequencies  and  the 
tension  are  slightly  modified  by  the  perturbations.  The  term  in  in  the 
tension  has  been  approximated  to  obtain  5.24b  and  c.  If  we  regard  X  and 
i>  as  Infinitesimal  and  omit  the  initial  conditions,  the  solution  of  5.24b 
and  5.24c,  with  only  the  driving  terms  kept,  is 


24  c 
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=  e|l  -  -^jjsinJl  5.25a 

=  5jlsin(e  +  g)  .  5.25b 


Substituting  for  X  in  the  last  term  on  the  left  of  5.24c, 
for  <p  on  the  right,  keeping  only  the  resonant  terms  yields  the  critical 
second-order  equation 


4(1  +  3j)n^X2  -  -18n^jelsin(2ll  +  g) 


5.26 


The  right  side  is  the  driving  resonant  term,  with  frequency  2(1  +  3/4j)n. 
This  differs  from  the  natural  frequency  on  the  left  by  a  term  of  order  j, 
so  the  factor  j  will  disappear  from  the  solution.  This  resonant 
displacement  has  the  value 


X^  =  -3elsln(2a  +  g) 


5.27 


If  e  =  .01,  I  =  1.5  deg,  the  amplitude  of  this  resonant  term  is  .045 
deg.  Hence,  although  the  resonance  exists,  its  effect  is  very  slight. 

The  same  type  of  analysis  can  be  applied  to  the  vibration  equations. 
However,  the  lack  of  a  driving  term  in  the  equation  analogous  to  5.24b, 
and  the  fact  that  the  lowest  vibration  frequency  is  thrice  rather  than 
twice  orbital,  allows  no  resonances  to  exist  to  at  least  second  order  in 
the  parameters.  Hence,  the  effect  of  oblateness  on  vibration  is  at  most 
to  change  the  frequencies  of  the  fundamental  modes  by  amounts  of  order  j , 
which  is  about  .00027  for  PACSAT  at  an  altitude  of  one  earth  radius.  We 
deduce  that  there  are  no  significant  effects  of  oblateness  on  vibration, 
and  the  only  effect  is  the  resonant  libration  5.27. 

We  next  consider  other  orbital  perturbations. 
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VI.  OTHER  ORBITAL  PERTURBATIONS 


In  this  section  we  shall  consider  the  effects  of  solar  and  lunar 
gravity  and  of  solar  radiation  pressure.  The  theory  has  been  developed 
by  previous  investigators.  Just  as  oblateness  causes  a  precession  of  the 
orbit  plane  around  the  equator,  solar  effects  produce  precession  around 
the  ecliptic,  and  lunar  gravity  produces  precession  around  the  lunar 
orbit  plane.  We  will  show  that  the  magnitude  of  the  precession  rates 
produced  by  the  sun  and  moon  is  very  small  compared  to  the  oblateness 
precession  rate,  and  infer  that  the  effects  on  llbration  and  vibration 
will  be  correspondingly  small. 

The  parameter  we  will  use  for  comparison  is  the  precession  rate  of 
the  node.  This  is  given  from  equation  5.7c  as 


3 

h  = 

R 


6.1 


For  PACSAT  at  an  altitude  of  one  earth  radius,  this  has  the  numerical 
value 


h  =  -.875  deg/day  . 


6.2 


The  corresponding  expressions  for  the  solar  and  lunar  precession 
rates  are [19] 


6.3 


where  D  denotes  the  disturbing  body,  of  mass  distance  R^,  and 

inclination  to  the  equator  Ijj.  For  the  sun,  we  have 
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3  X  10^ 


2  X  6.378  X  10“ 


1.495  X  10 


8 


=  8.53  X  10 


Ip  =  23.45  deg  . 


These  yield 


6 . 4a 


6.4b 


6 . 4c 


h  =  -^  yi  2.16  10^  X  3  X  10^  X  (8.53  x  10  x 

4 

=  -.00023  deg/day  . 


.762 


6.5 


Similarly,  for  the  moon,  we  have 


Mp  -  81.3 

E_  ^  2  X  6.378  X  10^ 
^  3.844  X  10^ 

18.3  <  Ip  <  28.5 


.0332 


6 . 6a 


6.6b 


6 . 6c 


where  the  third  of  equations  6.6  describes  the  precession  of  the  lunar 
orbit  plane  around  the  ecliptic.  Using  18.3  deg,  the  value  which 
produces  the  greatest  precession,  yields 
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2.16  X  lo’  ^ 


(.0332)  X  .852  =  -.00062  deg/day 


6.7 


Thus,  the  ratio  of  the  lunar  precession  rate  to  the  oblateness  precision 
rate  is  about  .0007,  and  we  infer  that  the  effects  on  libration  and 
vibration  will  be  reduced  by  about  the  same  relative  value.  Hence,  luni- 
solar  gravity  effects  will  be  negligible. 

We  next  consider  the  effect  of  solar  radiation  pressure.  This 
effect  has  caused  strange  behavior  on  the  part  of  past  satellites.  It 
produced  the  large  changes  observed  in  the  perigee  of  the  ECHO  balloon, 
and  an  interaction  between  solar  radiation  pressure  and  bending  produced 
by  solar  heating  caused  the  rapid  reduction  in  spin  rate  of  the  boom- 
carrying  satellites  Alouette  and  Explorer  XX[20] .  However,  we  shall  show 
that  for  PACSAT  the  solar  radiation  pressure  effects  are  negligible. 

The  solar  radiation  pressure  on  a  specularly  reflecting  sphere 
produces  the  acceleratlon[203 , 


JA-»- 

“cM° 


6.8 


where  J  is  the  solar  incident  flux,  A  the  projected  area,  c  the  velocity 
of  light,  M  the  mass,  and  a  the  unit  vector  directed  from  the  sphere 
toward  the  sun.  If  the  satellite  is  not  specularly  reflecting,  the 
magnitude  of  the  acceleration  is  reduced.  This  acceleration  is  in  the 
direction  opposite  to  the  solar  gravity,  and  may  be  thought  of  as  a 
reduction  in  the  effective  gravitational  constant.  For  PACSAT,  we  have 
the  numerical  values 


J  =  1.37  X  lo\/in^ 


6.9a 


c 


o 

3  X  10  m/sec 


6  •9b 


1 


6.9c 


A  ^  3 
M  4  rp 


.595  X  10  X  2.69  x  lO'" 


=  .0468  m  /kg 


from  which  we  deduce  the  acceleration: 


a  =  2.14  X 


10  ^ro/sec^ 


6.10 


The  acceleration  due  to  solar  gravity  is 


a  = 


6.11 


where  g  is  the  acceleration  of  gravity  at  the  earth's  surface.  This  has 
the  numerical  value 


a  =  9.8  X 


3  X  10^  X  (4.266  X  lo"^)^  =  5.35  x  lo‘^ 


6.12 


We  see  that  the  relative  reduction  in  the  effective  solar  gravity 
caused  by  solar  radiation  pressure  is  about  4  x  10  ^ .  Hence,  solar 
radiation  pressure  causes  a  very  small  change  in  an  already  small  effect 
and  can  safely  be  ignored. 

There  is  another  solar  radiation  pressure  effect  which  should  be 
considered.  The  satellite  will  pass  into  and  out  of  the  earth's  shadow. 
During  its  passage  into  the  penumbra,  the  only  effect  is  a  steady 
decrease  in  the  solar  flux.  The  decrease  is  uniform  along  the  satellite. 
However,  during  the  time  Interval  which  begins  when  the  lower  end  of  the 
satellite  enters  the  umbra  and  ends  when  the  upper  end  enters,  a  "wave  of 
darkening"  passes  up  the  satellite,  and  the  solar  radiation  pressure  is 
extinguished.  This  could  potentially  cause  both  libration  and  vibration. 
The  effect,  however,  is  negligible.  We  will  calculate  the  change  in 
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velocity  of  the  upper  end  of  the  satellite  relative  to  the  lower,  and 
show  that  this  change  is  far  below  detectability. 

We  will  calculate  the  effect  at  the  time  when  the  sun  is  on  the 
equator,  that  is,  the  equinox,  since  that  is  clearly  when  the  effect  is 
greatest.  The  configuration  is  depicted  in  Fig.  9.  The  sun,  of  diameter 
D  ,  is  at  distance  R  .  The  satellite  enters  the  penumbra  at  P,  the  umbra 

at  Q.  The  figure  is  not  to  scale,  and  the  lines  connecting  the  sun  to  P 

and  Q  should  be  nearly  vertical.  The  angle  If  which  measures  the  sun  s 
angular  diameter  is 

Y  =  ^  =  1.398  10  ^  ,00933  rad  =  .53  deg.  6.13 

\  1.495  X  lo“ 

For  PACSAT  at  one  earth  radius  altitude,  the  angle  0^  is  60  deg.  It 
is  easily  demonstrated  that  the  angle  POQ,  the  width  of  the  penumbra  as 
seen  from  the  center  of  the  earth,  is  also  equal  to  JT .  The  arc  length  PQ 

is  then  Ry  =  119  km.  The  velocity  of  the  satellite  in  orbit  is  5.59 

km/sec,  so  the  time  to  cross  the  penumbra  is  21.3  sec. 

If  we  draw  the  plane  from  the  satellite  instantaneous  position, 
measured  by  the  angle  8,  tangent  to  the  earth's  surface  and  extended  to 
the  sun,  and  neglect  atmospheric  refraction,  the  portion  of  the  sun's 
area  between  the  intersection  of  that  plane  and  the  left  edge  of  the  sun 
(Fig.  9)  is  observed.  It  is  easily  shown  by  the  law  of  cosines  and  a 
simple  calculation  that  the  solar  flux  is  reduced  from  its  "full  sun" 
value  J  to  the  effective  value 


J' 


T  2  .  ”1  r  YQ  /  ,1/2 

J  X  -Sin  [  (Sq  -  e)/Y  ] 


6.14 


At  the  time  when  the  lower  end  of  the  satellite  just  enters  the 
umbra,  the  upper  end  is  at  the  angular  position  0q  ■  L//3R.  The  flux  is 
reduced  to  the  amount : 
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Fig.  9  —  Passage  through  penumbra 
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J’ 


1.37  X  X  i 

7f 


74.4  w/m^ 


6.15 


and  the  acceleration  to  the  amount: 


a 


1 . 16  X 


10  m/sec 


6.16 


The  time  to  pass  into  the  full  umbra  is 


L 


X 


2ti 


.3  sec 


6.17 


where  T  is  the  satellite  period,  so  the  change  in  velocity  of  the  upper 
end  with  respect  to  the  lower  is 


Av  =  1.16 


X  10“®  X 


.3  =  3.5  X 


-9 

10  m/sec 


6.18 


This  is  far  below  any  detectable  amount,  and  also  much  below  any  error 
caused  by  emplacement,  so  the  effect  of  passing  into  or  out  of  the 
earth's  shadow,  and  thereby  losing  solar  radiation  pressure,  is 
completely  negligible. 
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VII.  MICROMETEOROID  IMPACT  EFFECTS 


All  the  external  forces  considered  to  this  point  have  very  large 
scale  factors.  The  gravity-gradient  and  oblateness  terms  have  scales  for 
variation  on  the  order  of  the  orbit  radius,  the  solar  and  lunar  effects 
scale  like  the  distance  to  the  sun  or  moon.  By  contrast,  micrometeoroid 
impacts  have  as  a  geometrical  scale  the  radius  of  the  impacting  body, 
which  is  small  compared  to  the  size  of  PACSAT.  Hence,  we  can  expect  the 
excitation  of  many  vibrational  modes  by  an  impact.  We  shall  assume  that 
the  impacting  body  is  sufficiently  small  that  it  does  not  sever  or 
produce  other  gross  effects  on  PACSAT. 

The  momentum  of  the  incoming  body  will  be  absorbed  in  a  region  of 
length  about  equal  to  the  meteoroid  diameter.  There  will  be  cratering 
and  spallation,  which  cause  ejection  of  mass  from  the  impacted  body.  The 
theory  of  this  effect  can  be  found  in  Ref.  21.  Approximately  half  the 
kinetic  energy  of  the  meteoroid  is  transferred  to  the  ejecta,  so  we  shall 
take  the  momentum  transfer  to  PACSAT  as  M^v^/Zz,  where  is  the  mass  of 
the  meteoroid  and  v^  is  its  velocity,  furthermore,  the  direction  of  the 
impacting  body  is  arbitrary.  The  component  of  velocity  transverse  to 
PACSAT  has  the  rms  value  so  the  transverse  component  of  momentum 

transfer  will  be  M^v^//6. 

The  theory  of  the  transverse  vibrations  of  struck  bars  was  developed 
by  Rayleigh  (Ref.  14,  p.  270),  and  the  theory  of  bars  in  orbit  is  almost 
the  same.  We  ignore  the  ellipticity  and  perturbing  forces,  and  let  the 
meteoroid,  with  effective  momentum  M^v^/Zd  and  diameter  d^,  strike  at  a 
point  Sq.  The  section  of  rod  about  of  length  d^  will  absorb  the 
momentum,  and  will  be  set  into  motion  with  an  initial  velocity  y^,  which 
is  determined  from  conservation  of  momentum: 


y^pdi  =  M^v^/i/6 


7.1b 
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Suppose  the  direction  o£  impact  is  in-plane.  The  theory  is  the  same 
for  out-of-plane  motion.  Let  the  initial  displacement  be  zero.  Then  the 
solution  of  the  vibration  equation  is  given  by  3.31,  where  =  0,  and 
Aj^(O)  is 


„  M.v. 


7.2a 


7,2b 


and  thus  the  complete  motion  is 


y(x,t) 


nL 


M.v, 


7.3 


where  the  libration  is  described  by  k  =  0,  Xq(x)  =  x, 
odd  vibratory  modes  are  to  be  included  in  7.3. 

tfe  calculate  the  motion  of  the  point  of  impact, 
over  possible  impact  locations,  we  obtain  the  simpler 


=  /S.  Even  and 

If  we  then  average 
expression: 


y(t)  = 


nL 


M.v. 


7.4 


The  series  which  appears  is  an  oscillatory  function  of  time,  whose  peak 
value  is  approximately  1.5.  Hence,  the  amplitude  of  the  transverse 


motion  is 
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Y  = 


M.v. 
1  1 

nLp 


7,5 


and  the  libration  is  about  3/8  of  this. 

Data  for  the  meteoroid  population  are  given  in  Ref.  21.  The 
incoming  velocity  is  taken  to  be  16  km/sec.  This  appears  rather  low, 
since  meteoroids  are  observed  with  incoming  velocities  up  to  70  km/sec. 

¥e  will  use  the  referenced  value,  but  observe  that  the  momentum  transfer 
scales  linearly  with  the  incoming  velocity.  With  the  appropriate 
constants  for  PACSAT,  the  vibration  and  libration  amplitudes  are  given  by 


Yy  =  6.05  X  loV(g)  -  7.6a 

Y  =  2.27  X  loSl.(g)  7.6b 

L  1 


in  centimeters.  If  we  had  taken  the  impact  at  the  tip,  these  values  are 
approximately  doubled. 

The  permissible  vibration  amplitude  is  .5  cm,  the  permissible  tip 
libration  amplitude  is  20  in  (1.5  deg).  The  impacting  mass  which  will 
produce  these  deflections  is 


=  8.3  X  10  ^g  (vibration)  7.7a 

=  .088g  (libration)  .  7  -  7b 


A  particle  of  mass  .088g  and  a  velocity  of  16  km/sec  has  a  kinetic  energy 
of  about  11,000  joules.  Such  an  impact  would  undoubtedly  smash  PACSAT  to 
the  proverbial  smithereens . 
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We  now  must  consider  the  probability  of  impact  of  particles  of  the 
indicated  size.  The  flux  of  meteoroid  particles  per  square  meter  per 
second  is  given  in  Table  I  of  Ref.  21,  which  we  reproduce  in  graphical 
form  in  Fig.  10.  We  see  that  the  flux  at  the  indicated  levels  is 


Flux 

(8.3  X  lO’^)  =  5.3  X  10‘® 

-2  -1 
m  s 

CO 

Flux 

(.088)  =  lO'^^ 

-2  -1 
ms  . 

7.8b 

Fig.  10  —  Meteoroid  flux  (from  Barengoltz) 
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The  projected  frontal  area  is 


A  =  N(iTr  +  )ld) 


7.9 


where  r  is  the  radius  of  the  sphere  (.595  cm),  Z  the  length  (.685  cm)  and 

d  the  diameter  (.132  cm)  of  the  rod,  and  N  is  the  number  of  sections. 

2 

This  works  out  to  a  frontal  area  of  9 . 62  m  .  The  product  of  this  and  the 
flux  gives  the  probability  of  impact  per  day  as 


Probability  (vibration)  =  .044  7.10a 

Probability  (libration)  =  .83  «  10  7.10b 


from  which  the  expected  time  between  impacts  is 


Interval  (vibration)  =  22  days  7.11a 

Interval  (libration)  =  1.2  x  10^  days  =  3300  years  .  7.11U 


If  the  representative  velocity  were  50  km/sec  instead  of  15,  the  interval 

for  vibratory  excitation  exceeding  the  permissible  would  be  only  8  days. 

Since  the  intervals  should  have  a  Poisson  distribution,  the  standard 

deviation  should  be  about  the  same  as  the  mean. 

We  should  check  the  consistency  of  our  assumptions.  If  the  incoming 

3 

particle  is  a  small  grain,  its  density  is  about  .6  g/cm  (Ref.  21).  The 

-6  3 

diameter  of  a  spherical  particle  of  mass  8.3  x  10  g  and  density  .6g/cm 
is  .025  cm,  small  compared  to  the  size  of  a  PACSAT  sphere.  Hence,  it  is 
reasonable  to  treat  the  momentum  absorption  as  concentrated.  The  kinetic 
energy  of  such  a  particle  is  about  1  joule,  and  the  initial  velocity  is 
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given  by  7.1b  as  about  2,S  m/sec.  The  PACSAT  structure  should  be  able  to 
sustain  such  an  impulse  without  rupturing. 

We  see  that  the  effect  of  raicrometeoroid  impact  on  PACSAT  vibrations 
is  very  Important.  Even  if  the  initial  alignment  is  perfect,  we  can 
expect  excessive  vibration  to  be  induced  by  meteoroid  impact  within  about 
a  20-day  interval.  This  is  much  too  frequent  for  satisfactory  operation. 
Some  method  of  damping  must  be  provided. 

Before  leaving  the  subject  of  meteoroid  impacts,  we  shall  consider 
the  excitation  of  longitudinal  vibrations.  Analogous  to  the  calculation 
of  the  transverse  momentum,  the  effective  longitudinal  component  of 
momentum  transfer  is  The  theory  of  longitudinal  vibrations  of 

bars  has  also  been  worked  out  by  Rayleigh  (Ref.  14,  p.  245).  The 
longitudinal  motion  z(s,t)  is  governed  by  the  differential  equation  and 
boundary  conditions: 


2  2 

d  z  ,  2  3  z 

— 2  =  b  — 2 

3t  3s 


7.12a 


at 


3z 

3s 


0 


b  = 


5.04  X  10  cm/sec  (steel) 


7.12b 


7.12c 


Here  o  is  the  volume  density  of  steel  and  E  is  Young's  modulus.  The 
parameter  b  is  the  velocity  of  sound  waves  in  steel.  We  choose  the 
initial  deflection  equal  to  zero,  and  the  initial  velocity,  analogous  to 
7.1b,  with  impact  at  s^. 


Zn  =  M.v.//5pd.  s.  -  .5d,  <  s  <  s„  +  .5d.  . 

Oil  'll  0  1  0  1 


7.13 
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The  solution  of  7.12a,  subject  to  the  boundary  conditions  7.12b  and 
initial  condition  7.13,  is 


z(s,t) 


M^v. 
1  1 


/6TTpb 


”  2Trks  2irksQ 
iCOS^; - cos- — 


sin27Tkbt/k 


7.14 


The  series  has  the  peak  value  if/2,  so  the  amplitude  for  longitudinal 
vibration  is 


M,v . 

1  1 

2/6Pt  ' 


7.15 


For  a  16  km/sec  particle,  this  becomes 


Z  =  .46M^(g)  cm  . 


7.16 


The  change  in  the  radar  cross  section  of  PACSAT  produced  by 
longitudinal  vibrations  is  not  given  by  the  analysis  of  Ref.  3  which  led 
to  Table  1,  and  we  will  calculate  it  here.  The  position  of  the  m  th 
sphere  is  given  by  tns  +  z(nis),  where  s  is  the  spacing  between  elements. 
Assuming  z  is  small,  if  we  define  the  mean  and  mean  square  values  of  z  by 


N/2 


z 


7.17a 


94 


N/2 


the  change  in  radar  cross  section  is 


Ai 

$ 


I2j^\2l2 


-2 

z 


)• 


—  -2 

The  moments  and  z  will  be  functions  of  time, 
is  at  the  center,  =  0,  then  the  moments  z  and 


M.v. 
r  1 

/grph 


k 


sin2Trkat 


2 

z  = 


Nz^ 


The  series  in  7.19a  has  the  amplitude  ti/2,  so  we 
cross  section 


A$/4>  =  (-^jyN . 


With  s  =  1.S75  cm, 


N  =  8  X  10  ,  and  z  given 


Al/f 


-5  2 

2.97  X  10 


7.17b 


7.18 


If  we  assume  the  impact 
_2 

z  are  calculated  to  be 


7.19a 

7.19b 

obtain  for  the  change  in 


7.20 

by  7.16,  this  gives 


7.21 
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Fot  a  1  percent  change  in  cross  section,  the  required  mass  is  18  grams. 

A  particle  sufficient  to  produce  a  transverse  vibration  of  .5  era  would 
produce  a  longitudinal  vibration  of  2.7  *  10  ^  cm,  clearly  a  negligible 
effect.  The  lowest  frequency  associated  with  longitudinal  vibrations  is 
b/L  =  3.3  sec  so  the  longitudinal  vibrations  are  very  small  and  very 
rapid  compared  to  the  transverse  vibrations. 

The  meteoroid  impact  can  also  produce  torsional  vibrations.  These 
will  not  have  any  effect  on  the  radar  cross  section,  so  we  need  not 
analyze  them.  The  effects  on  the  motion  should  be  comparable  to  the 
displacements  associated  with  longitudinal  vibrations. 

We  can  make  a  crude  estimate  of  the  force  involved  in  the  impact  of 
a  particle  which  produces  .5  cm  vibration  amplitude.  According  to  Ref. 
21,  a  meteoroid  impact  produces  a  hemispherical  crater,  with  a  radius 
given  by 


R 

c 


d. 

X 


1,  2,„.l/3 

2(aiV./S) 


7.22 


Here  a.,  v.,  and  d.  are  the  volume  density,  velocity,  and  diameter  of  the 

meteoroid,  and  S  is  the  yield  strength  of  the  material.  Assuming  the 

impact  is  on  a  sphere,  we  use  for  S  the  tensile  strength  of  aluminum, 

which  depends  on  the  forming,  annealing,  and  other  properties  of  the 

8  8  2 

spheres  and  should  be  between  ,4  x  10  and  2  x  10  Newtons/ra  .  With  the 

3  4  -4  X 

previous  numbers  =  600  kg/m  ,  v^  =  1.6  x  10  m/sec,  d^  —  2.5  x  10  m) , 

,7.22  yields  a  crater  radius  between  .1  cm  and  .17  cm.  We  use  . 14  cm  '  , 

as  representative.  The  crater  is  formed  by  a  fracturing  and  yielding 

process,  which  we  will  estimate  to  propagate  at  the  speed  of  sound  in 

bulk  aluminum  (6.4  x  10^  cm/sec).  Since  the  impact  is  hypersonic,  we 

should  actually  use  the  velocity  of  propagation  of  shock  waves,  but  this 

differs  from  the  sonic  velocity  by  a  factor  of  2  or  3  and  we  are  only 

interested  in  the  order  of  magnitude  of  the  force.  The  time  for  the 

meteoroid  to  be  brought  to  rest  we  estimate  by  dividing  the  crater  radius 

by  the  speed  of  sound,  and  obtain 
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At  =  R^/b  =  .14/6.4  X  10^  =  2.2  x  lo'^  sec  .  7.23 


The  average  force  is  the  momentum  change  divided  by  the  time 

Interval  At,  or 


F  =  M^v^//2At  =  8.3  X  lo'^  x  1.6  x  10^//2  x  2.2  x  lo'^  7.24a 


=  4  X  10  dynes  . 


7.24b 


We  see  that  this  force  is  very  large  compared  with  other  forces  that  may 
act,  but  it  has  a  very  short  duration. 

The  calculations  of  meteoroid  effects  have  been  based  on  the 
assumption  that  the  solid  sphere  stops  the  meteoroid.  If  PACSAT  were 
built  of  hollow  spheres,  the  meteoroid  might  punch  its  way  through  and 
exit  carrying  off  some  of  the  momentum.  We  take  as  a  reasonable 
approximation  that  the  meteoroid  should  have  sufficient  velocity  to  punch 
two  craters  of  depth,  given  by  7.22,  equal  to  the  thickness  of  the  sphere 
wall.  The  tensile  strength  of  plastic  should  be  used,  which  we  expect  to 
cover  the  same  range  as  does  aluminum  (phenolic  plastics,  the  most  up-to- 

g 

date  type  in  the  critical  tables,  cover  the  range  .7  to  1.7  x  10 
2 

Newtons/m  ).  According  to  the  immediately  preceding  calculation,  an  8.3 
X  lo'^g  particle  will  just  penetrate  a  shell  with  wall  thickness  .1  cm, 
and  will  certainly  not  have  enough  energy  remaining  to  penetrate  the 
shell  a  second  time.  Hence,  if  the  hollow  spherical  shell  has  a 
thickness  of  one  millimeter,  which  is  certainly  reasonable,  the  analysis 
for  the  solid  sphere  should  still  be  correct. 

We  conclude  that  meteoroid  Impacts  can  have  a  very  important  effect 
on  the  usage  of  PACSAT,  since  impacts  sufficient  to  cause  a  radar  cross- 
section  loss  of  1  dB  can  be  expected  to  occur  about  every  20  days.  A 
loss  of  10  dB  can  be  expected  about  every  50  days.  The  oscillations 
associated  with  these  impacts  must  be  damped  if  PACSAT  is  to  operate 
satisfactorily. 

The  next  and  final  topic  for  investigation  is  the  effect  of  bending 
caused  by  thermal  stresses. 
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VIII.  THERMAL  BENDING  EFEECTS 


There  have  been  numerous  investigations  of  the  bending  of  long  booms 
under  solar  input.  Most  of  these  studies  have  been  applied  to  the  Radio 
Astronomy  Explorer  satellite  (1968-55A),  which  carries  four  750-foot 
booms  in  a  cruciform  arrangement  around  the  satellite  body,  plus  a  damper 
boom.  The  booms  are  extendible  beryllium-copper  thin-walled  tubing  with 
overlapping  portions.  A  study  by  Yu{22],  with  comments  by  Augusti[23], 
Jordan[24],  and  Yu[25j,  concerning  the  stability  of  a  boom  making  a  fixed 
angle  with  the  sun- line,  showed  after  much  disagreement  that  a  boom  is 
stable  when  the  angle  between  it  and  the  sun- line  is  less  than  90  deg, 
but  is  unstable  if  the  angle  exceeds  90  deg.  Sikka  et  al.[261  determined 
the  temperature  distribution  and  curvature  in  long  solid  cylinders  in 
space,  tfe  shall  use  their  results  directly.  Thermal  curvatures  in  long 
booms  have  been  studied  by  Eby  and  Karara[27],  Vigneron[ 28] ,  and 
Graham[29],  with  particular  attention  paid  to  the  effects  of  thermal  time 
constants  on  the  motion. 

The  problem  to  be  treated  may  be  visualized  as  follows.  PACSAT,  in 
orbit  around  the  earth  in  a  nominally  vertical  position,  receives  direct 
heat  from  the  sun  and  reflected  heat  from  the  earth.  The  reflected  heat 
causes  temperature  variations  along  the  rod  which  would  primarily  lead  to 
longitudinal  displacements  and  vibrations,  which  we  may  neglect  per  the 
investigation  on  micrometeoroids.  The  direct  heat  causes  temperature 
variations  across  the  rod.  The  rod  expands  more  on  the  heated  side  than 
on  the  unheated,  which  causes  it  to  curve  away  from  the  sun-line.  This 
thermal  curvature  will  couple  to  the  mechanical  curvature  used  in  all 
preceding  sections  of  this  report,  and  we  can  expect  substantial 
excitation  of  vibratory  modes.  In  particular,  if  the  solar  input 
contains  third  harmonics  of  the  orbital  frequency,  there  will  be  resonant 
excitation  of  the  lowest  vibrational  mode,  and  large  deflections  may 
result.  This  proves  to  be  the  case,  and  thermal  bending  effects  will 
apparently  render  the  present  design  of  PACSAT  unfeasible  for 
satisfactory  communications. 
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We  begin  by  summarizing  the  investigation  of  Sikka  et  al.[26].  They 
treat  a  long  solid  cylinder  of  homogeneous  isotropic  material.  The  flux 
from  the  sun  is  incident  normally,  and  a  fraction  a  is  absorbed.  Thermal 
variations  along  the  cylinder  are  neglected.  The  steady-state 
temperature  is  found.  This  temperature  distribution  satisfies  Laplace's 
equation,  with  the  boundary  condition  which  matches  the  thermal  input 
with  the  losses  by  conduction  to  the  interior  and  radiation  to  space.  In 
terms  of  the  radial  coordinate  r  from  the  center  of  the  cylinder  and  the 
angular  coordinate  ^5  around  the  cylinder,  measured  from  the  sun-line,  the 
temperature  distribution  is  approximately  given  by 


1^ 

/I  a  J \4  ,  1  aJ  ^  _  _ 

where  a  is  the  absorptivity,  £  the  eraissivity,  J  the  solar  input  [1.37 
X  lo\/m^],  0  is  Stefan's  constant  [5.67  x  10  ®w/m^/ (“K)^] .  and  k  is  the 
thermal  conductivity.  Terms  of  higher  order  angular  dependence  do  not 
contribute  to  the  curvature.  For  a  gray  body  whose  absorptivity  and 
emlssivity  are  equal,  the  first  term  in  8.1  is  296°K.  For  a  steel  rod, 
k  =  16  w/m/^K,  r  =  .066  cm,  the  second  term  is  .O28acos0,  so  even  for  a 
perfectly  absorbing  rod  the  variation  in  temperature  around  the  rod  is 
much  less  than  0.1*^K. 

The  solution  of  the  thermoelastic  problem  of  the  bending  of  a 
cylinder  under  thermal  expansion  is  well  known[30].  The  curvature  K 
produced  in  the  cylinder  due  to  the  temperature  distribution  T  is  the 
same  as  the  curvature  induced  in  a  long  cylinder  due  to  application  of 
stress  eET  at  the  ends,  where  e  is  the  coefficient  of  thermal  expansion 
and  E  is  Young's  modulus.  When  the  stress  is  calculated  with  the 
temperature  distribution  1.1,  the  simple  result  from  Ref,  26  for  the 
curvature  is 


2. 

"2k 


8.2 
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For  a  steel  rod,  the  coefficient  of  expansion  e  is  1.8  x  10  /  K,  so  the 

curvature  is 


K  =  7.7  X  10 


-4 

a 


-1 

m 


8.3 


and  the  radius  of  curvature  Rj,,  the  reciprocal  of  the  curvature,  is 


=  1.30/a  km  . 

Xv 


8.4 


It  is  apparent  at  this  stage  of  the  analysis  that  PACSAT  is  in. 
trouble.  The  requirements  on  straightness  indicate  a  radius  of  curvature 
of  50,000  km,  which  would  necessitate  an  absorptivity  of  2.6  x  10 
(.0026  percent)  if  the  thermal  radius  of  curvature  8.4  represents  the 
actual  curvature.  This  absorptivity  appears  to  be  beyond  the  current 
state  of  the  art.  We  shall  find  the  actual  radius  of  curvature,  and  will 
show  that  it  leads  to  an  absorptivity  requirement  of  .07  percent,  which 
is  also  beyond  the  state  of  the  art. 

Equation  8.2  gives  the  curvature  when  the  rod  is  straight  and  the 
sunlight  is  perpendicular  to  the  axis.  When  the  rod  is  bent,  and  the  sun 
is  in  some  other  direction,  the  curvature  should  be  perpendicular  to  the 
slope,  be  in  the  plane  containing  the  slope  and  the  sun- line  vector,  be 
directed  away  from  the  sun,  and  be  proportional  to  the  cosine  of  the 
angle  between  the  normal  and  the  sun-line.  These  conditions  are 
satisfied  by  the  expression[6] 


8.5 


where  K  is  expression  8.2,  3r/9s  is  the  slope,  and  o  is  the  unit  sun¬ 
line  vector.  We  have  used  the  unit  magnitude  of  3r/3s. 
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The  shadows  cast  by  the  spheres  will  shield  portions  of  the  rods 
from  solar  input.  The  configuration  is  shown  in  Fig.  11.  The  satellite 
is  in  equatorial  circular  orbit  and  is  approximately  radial.  In 
equatorial  coordinates,  the  slope  is  approximately 

—  ej  ^  =  i  COS0  +  1  sin0  8.6 

ds  R  X  y 


Fig.  11  —  Shadowing  of  rods  by  spheres 
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while  the  sun- line  vector  a  is 


0  =  i  COS0  +  i  cosesinO  +  i  sinEsin0  ' 

X  s  y  s  z  s 

where  £  is  the  inclination  to  the  ecliptic  (23.45  deg)  and  0^  is  the 
solar  longitude  measured  from  Aries.  Ellipticity  of  the  earth  s  orbit 
around  the  sun  has  been  neglected..  The  angle  3  of  Fig.  11  is  given  by 


CDS?  =  —  •  CT  =  CDS0  COS0  +  cosesinB  sln9  .  8.8 
3s  s  s 

With  the  length  of  the  rod  equal  to  £,  and  the  radius  of  the  sphere  equal 
to  r,  it  is  easily  seen  that  the  length  of  the  unshadowed  part  of  the  rod 
is  given  by 


Z  =  I  -  r(csc3  ■  1)  • 
u 


The  radius  of  the  rod  has  been  ignored,  so  the  shadow  is  regarded  as 

going  straight  across  the  rod.  For  PACSAT,  Z  ~  .685  cm,  r  =  .595  cm,  and 

for  Z  to  be  positive  3  must  be  greater  than  27.7  deg.  Thus,  during  the 
u 

portions  of  the  orbit  when  PACSAT  is  sufficiently  close  to  being  directly 

between  the-  earth  and  the  sun,  or  for  those  when  PACSAT  is  nearly, 

opposite  the  sun,  there  will  be  no  thermal  bending  effects.  The 
opposition  condition  will  put  PACSAT  into  eclipse.  If  the  altitude 
exceeds  1.15  earth  radii,  PACSAT  will  be  in  full  internal  shadow  before 
it  goes  into  eclipse.  We  shall  consider  this  condition  only. 

'pjjg  thermal  curvature  K  is  to  be  added  to  the  mechanical  curvature 
to  form  the  strain  energy.  Thus,  we  have  the  new  form  for  the  strain 
energy  2.16; 
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V 


B 


rEI/ds 


3^r 


as 


K-^x 

as 


8.9 


The  vector  ^  has  a  maximum  variation  L/Rq,  the  ratio  of  the  satellite 

^  “  8 

length  to  the  distance  to  the  sun.  This  is  numerically  10  ,  so  0  will 

be  treated  as  independent  of  s . 

When  the  expression  in  brackets  in  8.9  is  squared  (squared  means  the 
square  of  the  magnitude),  the  square  of  the  first  term  gives  the  previous 
results.  The  square  of  the  second  term  may  be  shown  to  be  equivalent  to 
a  10~®  relative  change  in  the  tension.  For  the  cross-term,  we  may 
simplify,  since  the  curvature  is  perpendicular  to  the  slope,  and  obtain 


-EIK/ds 


a^r 

as^ 


a 


8,10 


forming  the  variation  and  integrating  by  parts, 


2 

-EIK/dsCT  •  — ^6r 

as 


-EIK  o  .  -Sr 


L 

2 

_L 

2 


rj 

fds-— 

ds 


8.11a 


8.11b 


The  second  term  is  zero,  since  a  is  independent  of  s.  Hence,  the  thermal 
curvature  does  not  affect  the  differential  equation  2.19.  The  first  term 
in  8.11b  must  be  combined  with  the.  corresponding  term  in  2.18. 
Consequently,  we  obtain  the  new  boundary  conditions: 
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At  s  =  ±7 


— y  =  Ka 


3s^ 


=  0  . 


8.12a 


8.12b 


These  equations  still  correspond  to  no  bending  or  shearing  at  the  ends. 
We  see  that  there  are  no  thermal  bending  effects  on  the  libration.  This 
is  obvious  anyway,  since  an  induced  curvature  cannot  to  first  order 
affect  a  rigid-body  motion. 

Before  proceeding  to  solve  the  vibration  equations,  we  need  an 
estimate  of  the  validity  of  the  use  of  the  steady-state  thermal  theory. 
The  thermal  time  constant  of  a  sphere  or  rod  Is  given  by 

2 

c  or 


where  c  is  the  specific  heat,  0  the  volume  density,  r  the  radius  of  the 
1  ^ 

sphere  or  rod,  and  k  the  thermal  conductivity.  This  expression  gives  for 
the  aluminum  sphere  and  the  steel  rod  the  values  .3  sec  and  .1  sec.  The 
vibratory  motions  and  the  rate  of  variation  of  the  solar  input  both  have 
time  scales  of  many  minutes,  so  the  condition  of  steady-state  temperature 
is  well  satisfied. 

Assume  a  circular  equatorial  orbit  and  form  the  equations  of 
vibration.  The  complete  rod  participates  In  the  bending,  but  only  the 
unshadowed  part  experiences  thermal  input.  Thus,  when  we  introduce 
normalized  variables,  a  factor  involving  should  appear  in  the  boundary 
conditions,  but  not  in  the  differential  equation.  The  result  for  the 
iu-plane  vibration,  analogous  to  3.27,  is 


3^2  = 


3x 


2  ®”2  ’"'“2 


8.14a 
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at  X  =  ±1 


=  L  Ka^/A  =  Q 


0 


8.14b 


8.14c 


where  component  of  the  sun-line  vector  in  the  in-plane 

direction.  We  obtain  readily 

0  =  -COS0  sin9  +  costsinB  cos0  . 

2  s  s 


Noting  the  definition  of  Q  in  8.14b,  we  solve  the  system  8.14  by  the 
trial  form 


=  "(3x^  -  1)Q  +  V  . 


8.16 


The  first  factor  has  been  chosen  to  satisfy  the  boundary  conditions  8.14b 
and  8.14c  and  to  have  vanishing  average  value.  Hence,  V  will  satisfy  the 
homogeneous  form  of  the  boundary  conditions  8.14.  The  differential 
equation  for  V  becomes 


—V"  =  —(1 

3  ax'- 


2,3V 
X  )— 


4 

ax 


-  (3x^  -  i)(|q"  +  q)  • 


3.17 


The  solution  may  be  expanded  in  the  form 


V 


(e)Xj^(x) 


8,18 
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which  satisfies  the  boundary  conditions,  and  leads  to  the  differential 
equation  for  Vj^: 


''1 


“^k 


=  --^(Q"  +  9Q) 


£ 


(3x^ 


-  l)Xj^(x)dx 


8.19 


Let  the  effective  solar  input  be  expanded  into  the  as  yet  unspecified  set 
of  frequencies  \>^  by 


Q  =  ZQ^sin(v^0  +  . 


8.20 


Then  the  solution  of  8.19  is 


\  = 


(Sx'^  -  l)X^(x)dx 


9  - 

- ^sin{v^6  +  (fi^)  8.21a 

k  1 


(3x' 


l)Xj^  (x)dx 


9  -  w" 

Q  +  ■  -  ^sin(v^e  +  (fi^) 

U).  -  V  . 

k  1 


.  8.21b 


When  this  solution  is  substituted  into  8.18,  and  the  result  into 
8,16,  the  term  in  Q  in  the  second  bracket  in  8.21b  leads  to  the 
cancellation  of  the  first  term  in  8.16.  The  -1  in  the  integral  does  not 
contribute,  since  all  Xj^  have  vanishing  average  value.  The  numerical 
values  of  the  Integrals  are  found  to  be 

8.22a 

8.22b 


(x)dx  =  .4216 


x)dx  <  .0005  k  >  1  . 
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Hence,  the  sum  over  k  in  8.18  reduces  to  a  single  term.  We  have 


2 

“l 

=  .2108X^(x)5:q^~ 


9 


— ssln(v,6  + 
2  1 
V , 

1 


8.23 


2 

Since  to^  =  9.00645,  we  see  that  the  deflection  is  small  unless  is 
near  3.  Hence,  the  input  Q,  which  depends  on  the  orbital  and  solar 
frequencies,  should  be  expanded  in  Fourier  series,  and  those  terms  near 
thrice  orbital  frequency  retained.  Thus,  we  represent  the  resonant  part 
of  Q  in  the  form 


03 

Q_  =  FQ.sin(3S  +  i9  +<(),)  . 
R  r  s  1 

— oo 


8.24 


The  calculation  of  the  coefficients  is  very  laborious.  They  are 

represented  as  double  Integrals,  which  must  be  expanded  as  slowly 

.2  .2 

convergent  series  in  the  ratio  sin  e/sin  S  =  .733.  Rather  than 
“  max 

present  the  details,  we  simply  give  the  result.  Using  the  numerical 
value  of  K  from  8.3  and  L  =  548  m,  the  deflection  is 


w 


2 


aX^(x) 


sln(30  -  39^)  (2.5078  +  .6058cos2e^ 


+  cos(3G  -  36^) (  .55O9sin20^ 


+  . O374cos40  ) 

s 

+  .O469sin40^) 


8.25 


As  mentioned  at  the  end  of  Section  III,  the  "throw’  of  that 

is,  the  difference  between  its  maximum  and  minimum  values,  is  2.341. 
This  enables  us  to  define  the  deflection  amplitude  as 
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W  =  aF(0g)sin[30  -  38^  +  §(0^)] 


8.26 


where  F  and  f  are  the  amplitude  and  phase  of  2.341  times  the  bracketed 
expression  in  8.25.  We  plot  F(0^)  and  ^(0^)  in  Fig-  12.  They  have  a 
period  of  one-half  year.  As  can  be  seen  from  the  figure,  §  differs  from 
zero  by  only  about  14  deg,  but  F  varies  from  about  7.5  m  at  the  equinox 
to  4.5  m  at  the  solstice. 

We  see  that  if  W  is  to  be  less  than  .5  cm,  a  must  be  less  than 
.0007.  The  excitation  of  the  resonant  mode  is  caused  by  shadowing,  which 
induces  the  third  harmonic  in  o^.  If  there  were  no  shadowing,  the  third 
harmonic  would  be  produced  by  eclipsing.  Incidentally,  there  is  no 
resonant  excitation  in  the  out-of-plane  motion. 


Time  (days) 

Fig.  12  —  Annual  variation  of  thernnal  deflection 
W  =  aF  (ej  sin  [3  [6 -0^) +<i>{  9,^] 
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The  requirement  that  a  be  less  than  .0007  is  probably  impossible 
with  the  current  state  of  plating  and  coating  materials.  Silver  plating 
achieves  an  a  of  only  .1  (Refs.  27,  29).  A  fresh  mirror  surface, 
produced  by  evaporating  silver  on  glass,  has  absorptivity  of  .01  -  .02 
over  the  solar  spectruin[31] .  It  is  highly  unlikely  that  this  surface 
could  be  formed  on  a  metal  substrate,  even  with  heroic  polishing. 

It  is  possible  to  improve  the  situation  somewhat  by  making  the  rod 
of  Invar  (34  percent  Ni  steel),  the  standard  low-expansion  material. 
Unfortunately,  the  curvature  depends  on  the  quotient  of  the  expansion 
coefficient  and  the  thermal  conductivity,  and  these  tend  to  vary  in  the 
same  manner.  For  steel,  e  —  16  x  10  ,  k  =  16,  and  the  quotient  is 

1.1  X  10  For  Invar,  e=1.5x  10^,  k=ll,  and  the  quotient  is 
1.4  X  10  an  Improvement  of  about  a  factor  of  9.  This  leaves  a  required 
ot  of  .006,  still  beyond  the  state  of  the  art,  and  says  nothing  about  the 
difficulty  of  working  Invar. 

The  use  of  plastic  or  ceramic  materials  would  make  matters  worse. 
Their  expansion  coefficients  are  somewhat  lower  than  metals,  but  their 
thermal  conductivities  are  much  lower.  Data[32],  with  the  necessary 
changes  of  units,  show  that  for  vinyl  acetate  rigid  copolymer, 
e  =  7  -  18  X  10  k  =  .08,  e/k  =  88  -  226  x  10  for  phenolic  plastics, 
e  =  30  -  40  X  10‘^,  k  =  .09  -  .16,  e/k  =  190  -  440  x  10~^,  and  for  glass 
e  =  .8  -  14  X  lO'^,  k  =  .05,  e/k  =  16  -  280  x  lo'^.  All  of  these  are 
considerably  worse  than  steel.  It  appears  that  little  can  be  done  by 
choice  of  materials  to  improve  the  thermal  bending  situation. 

It  appears  from  equation  8.26  and  the  numerical  data  that,  under  the 
best  of  circumstances,  thermal  bending  will  cause  PACSAT  to  oscillate 
with  a  period  of  about  80  minutes  and  an  amplitude  of  several 
centimeters,  far  too  great  for  satisfactory  communications.  There  may  be 
some  possible  redesigns.  If  a  thermal  shield  could  be  put  around  the 
rods,  it  might  be  possible  to  apply  coatings  in  cascade  to  reduce  the 
absorptivity.  Also,  if  the  satellite  is  built  with  flexible  sections,  so 
that  the  curvature  is  not  maintained  along  the  length,  it  would  present  a 
scalloped  appearance  and  have  a  much  smaller  amplitude  of  oscillation. 

As  was  shown  by  Burke{ll],  the  "sticky  hinges"  proposed  by  Stanford  will 
not  work  for  this  purpose,  since  they  will  lock  and  allow  the  curvature 
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to  propagate.  The  flexible  portions  must  transmit  the  gravity-gradient 
tension,  or  the  whole  stabilization  system  will  not  function.  The 
engineering  of  such  concepts  is  beyond  the  scope  of  the  present 
investigation  (a  euphemism  for  saying,  "l  don't  know  how”). 
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IX.  CONCLUSIONS 


Although  the  analysis  in  this  report  has  been  quite  lengthy,  since 
there  were  many  topics  to  cover,  the  conclusions  are  quite  brief: 

1.  The  verticality  and  straightness  conditions  must  be  satisfied 
initially,  which  will  be  at  least  an  extremely  difficult 
deployment  task. 

2.  The  major  causes  of  flexure  are  micrometeoroid  impacts  and 
thermal  bending.  Micrometeoroid  impact  may  produce  excessive 
vibration  at  tens  of  days  occurrence  times.  Thermal  bending 
effects  will  excite  resonant  oscillations  of  large  amplitude  and 
are  so  severe  that  it  is  highly  unlikely  that  they  can  be 
overcome . 

3.  The  flexural  misbehavior  of  PACSAT  in  orbit  is  such  that  it  is 
most  improbable  that  the  existing  design  of  PACSAT  can  perform 
its  communications  function. 


Ill 
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